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DEFINITION OF SYMBOLS 


Symbol 

Definition 

Units 

A 

Multiple of mass considered, 

( - ) 


Area 

in2 

C 

Specific heat 

cal 

mole°K 

E 

Energy 

cal 

g 

F 

Any function 

( - ) 

G 

Gibb's function 

cal 

H 

Enthalpy 

cal 

I 

Impulse 

sec 

K 

Equilibrium constant 

( - ) 

L 

Atomic weight 

{ - ) 

M 

Molecular weight 

g 

mole 

Ma 

Mach number 

{ - ) 

N 

Multiple of mole numbers (A • n = N 

mole 

P 

Pressure 

atm 

Q 

Heat energy 

cal 

R 

Specific gas constant 

cal 

R 

Universal gas constant, 1.98726 

cal 

mole°K 

S 

Entropy 

cal 

°K 
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DEFINITION OF SYMBOLS (Cont'd) 


Symbol 

T 

U 

V 
W 
X 

Y 
Z 
a 

b 

c 


d 

e 

g 

go 

h 


Definition 

T emperature 

Internal energy 

Volume 

Work energy 

Element 

Element 

Element 

Number of gram atoms of element Z 

Number of gram atoms of element Y 

Number of gram atoms of element X 

Concentration 

Specific heat per unit mass 

Differential 

Euler number, Z. 7182818 
Gibb's function per unit mass 

Gravitational constant, 32.1740 

Enthalpy per unit mass 


Units 
°K 
cal 
ft^ 
cal 
( - ) 

( - ) 

( - ) 
g ■ 

mole 

g 

mole 

g 

mole 

mole 

ft' 

cal 

g°K 

( - ) 

( - ) 

cal 

g 

ft 

sec' 

cal 

g 
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DEFINITION OF SYMBOLS (Cont*d) 


Symbol 

m 

n 

q 

s 

u 

V 

w 

X 

y 

z 

a 

y 

6 

A 

a 

e 


Definition 

Units 

Mass, slope of a curve 

g 

Mole number 

mole 

Heat energy per unit mass 

cal 

g 

Entropy per unit mass 

cal 

g°K 

Internal energy per unit mass 

cal 

g 

Specific volume 

fti 

lb 

Flow velocity 

ft 

sec 

Variable 

( - ) 

Variable 

( - ) 

Geodetic height 

ft 

Variable 

( - ) 

Angle 

o 

Specific heat ratio 

( - ) 

Difference 

( - ) 

Difference 

{ - ) 

Partial differential 

( - ) 

Expansion ratio 

( - ) 

Number of oxidizer atoms 

( - ) 
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DIFINITION OF SYMBOLS (Cont'd) 


Symbol 

Definition 

Units 


Mixture ratio 

( - ) 

V 

Valence 

( - ) 

1 

p 

Number of fuel atoms 
Density 

( - ) 

lb sec 2 
ft-* 


Stochiometric mixture ratio 

( - ) 


SUBSCRIPTS 


a 

Number of gram atoms of element Z 

mole 

b 

Number of gram atoms of element Y 

g 

mole 

c 

Number of gram atoms of element X 

g 

mole 

ch 

Chamber 

( - ) 

Des 

Desired 

( - ) 

equ 

Equilibrium 

( - ) 

F 

Fuel 

( - ) 

f 

Formation 

( - ) 

g 

Gas 

( - ) 

i 

Pertaining to the i^^ element or compound 

( - ) 

K 

Kinetic 

( - ) 

1 

Liquid 

( - ) 


X 



DEFINITION OF SYMBOLS (Cont'd) 


Symbol 

Definition 

Units 

max 

Maximum 

( - ) 

O 

Oxidizer 

( - ) 

P 

Potential, at constant pressure 

( - ) 

ref 

Reference 

( ) 

s , 

At constant entropy, sound 

( - ) 

sp 

Specific 

( - ) 

st 

Stagnation 

( - ) 

stat 

Static 

{ - ) 

T 

At temperature T 

( - ) 

f 

Throat 

{ - ) 

tot 

Total 

{ - ) 

V 

At constant volume 

{ - ) 

vap 

Vapor, vaporizing process 

( - ) 

X 

Element 

( - ) 

Y 

Element 

( - ) 

Z 

Element 

( - ) 

o 

Correct solution 

( - ) 

03 

Any location in the flow 

( - ) 



Symbol 

o 
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DEFINITION OF SYMBOLS (Cont'd) 
SUPERSCRIPTS 
Definition 

At reference level of one atmosphere 

Molar quantity 

Derivatives 


Units 
( - ) 
( - ) 
( - ) 
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CALCULATION OF ROCKET PERFORMANCE PARAMETERS 
BASED ON THE EQUILIBRIUM COMPOSITION OF 
THE COMBUSTION PRODUCTS 


By Klaus W. Gross 

George C. Marshall Space Flight Center 


SUMMARY 

The theory and development of equations are presented for calcu- 
lating the composition of combustion products for frozen and shifting 
equilibrium, and a combination of both, based on a chosen temperature 
limit. The equilibrium compositions of combustion products are used to 
calculate rocket performance parameters based on a one-dimensional flow 
model. 

Furthermore, a technique is described for determining the maximum 
specific impulse for a constant chamber pressure. A procedure is also 
included for establishing the optimum chamber pressure and propellant 
mixture ratio for use when the thrust level of a specific rocket engine 
is uprated. 


INTRODUCTION 

Theoretical equations are developed for the solution of equilibrium 
compositions of combustion products obtained in liquid propellant 
rocket motors by minimizing the Gibbs function. Subsequently, the more 
important rocket performance parameters can be calculated, such as: 
characteristic velocity, specific impulse, thrust coefficient, tempera- 
ture, expansion ratio, thrust, flow rate, stochiometric mixture ratio, 
molecular weight of the combustion products, isentropic exponent, local 
velocity, local speed of sound, specific heats at constant pressure and 
volume, specific gas constant of the combustion products, and local 
static pressure. The indicated performance parameters can be deter- 
mined for desired local static pressures and expansion ratios. 

A maximum specific impulse calculation routine for constant chamber 
pressure is presented as well as an optimum chamber pressure and mixture 
ratio calculation to obtain a desired thrust level. The thermodynamic 
properties, entropy, enthalpy and specific heat at a pressure of one 
atmosphere are determined from polynomials for a temperature range from 
298.15°K to 5000°K. 


The calculation procedure is limited as to number of elements 
taking part in the combustion and formed combustion compounds. However, 
space is provided to include or exchange additional elements and 
compounds. The basic program uses only a portion of available computer 
storage space. This allows the combination of this procedure with 
other programs. Many programs have been written that solve for 
equilibrium compositions automatically. This program uses mainly the 
approach of Gordan et al. (Ref. 1). However, changes have been made, 
and additional computation procedures are provided. 

Grateful acknowledgement is given Mr. Richard Lewis of General 
Electric Company assigned to MSFC Computation Laboratory. Without his 
proficiency in mathematics and machine techniques this effort might 
well have been far less productive. 
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CONSERVATION OF ENERGY 


Energy exists in various forms, e.g., kinetic energy, potential 
energy, electrical energy, chemical energy, heat energy, etc. Experi- 
ence has shown that in a closed system the sum of all energies does not 
change with time. An addition of heat energy (AQ) to a system can cause 
change in all existing energies. 


AQ = AU + AW + AEp + AEk + 2AE (1) 

The energy differences between state 1 and 2 per unit mass can be 
expressed as follows: 


1. Internal energy 

2. Work energy 


AU 

na 

AW 

m 


= U2 - ui 


= P2V2 - PiVi 


P is the pressure and v the specific volume 

AEp 


3. Potential energy 


m 


= Z2 - zi 


z is the geodetic height 

AEk ^ 2 ^ ^1^ 

4. Kinetic energy = 

w is the flow velocity and gQ is the gravitational constant. 
Neglecting the other forms of energy, the first law of thermodynamics 


reads 

W2^ - 

q 2 - qi = U2 - Ux + — + 


Z2 


Zj + P2 V2 - Pi Vi 


( 2 ) 
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CONSTANT ENTHALPY PROCESS 

The enthalpy per unit mass(h) is defined as the sum of the internal 
energy per unit mass (u) , plus the product P • v, the work energy required 
to transport the masses. This amount of work, however, no longer exists 
in the considered mass, but has been transferred to the surroundings. 

With h = u + Pv (3) 

qz - qi = h2 - hx + 3^- +22-^1 

The combustion process in a rocket combustion chamber can be assumed 
to be adiabatic, which means no heat will be added or taken out of the 
system and - q^= 0. The injection velocity of the propellants and 
also the velocities of the formed combustion products in the combustion 
chamber are relatively small and can be neglected. Therefore, 

^ 

2go 

Futhermore, the potential energy difference can be neglected. For a 
comubstion chamber with a flow in the vertical direction, the difference 
in height between the inlet and outlet cross sectional areas is small, 
and for a horizontal chamber the height differences between the centers 
of geometry of the two areas is zero. Therefore, 

Z 2 - zi = 0 

This leaves only the two enthalpy terms before and after the combustion, 
and they are equal, 

hi = h2 (5) 

Ui + Pi Vi = U2 + P 2 V2 
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ENTROPY 


The entropy change of an ideal gas can be developed from the 
second law of thermodynamics , 

Tds = dh - vdP (6) 

The enthalpy change (dh) can also be expressed by 
dh = cp dT 

Applying the ideal gas equation 

Pv = RT (7) 




we obtain 


ds 



( 8 ) 


Considering an isothermal process with dT = 0 and integrating the 
equation yields 


S 2 - Si = -R(lnP 2 - In Pi) (9) 

The third law of thermodynamics states: The entropy of matter is 

equal to zero tor O^K. Therefore, the entropy represents absolute values 
when this point is used as a reference level. Furthermore, if a pressure 
reference of one atmosphere is selected, we obtain 
S 2 - Si = Rln P 2 

Entropy of a gas per unit mass can now be calculated for a certain 
temperature and pressure 

- RlnP (10) 
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The Gibbs - Dalton law indicates that the entropy of a nonreactive 
mixture of ideal gases is equal to the sum of the entropy values of the 
individual constituents, considering their quantities present. Since 
most thermodynamic tables list the molar entropy of a species, the total 
entropy of a mixture can be calculated when the individual mole numbers 
of the compounds are known: 

s = S ni(s^)i (10 

i 

CONSTANT ENTROPY PROCESS 

The expansion of the combustion products through a De Laval 
nozzle occurs very rapidly and almost without friction losses or heat 
exchange with the surroundings. This is nearly equivalent to a reversible 
adiabatic process. Therefore, the total entropy of the gas mixture will 
hardly change during the expansion process, and the assumption of a 
constant entropy process is justified. 


s = ^ch 


( 12 ) 


Knowing the static pressure at a specific location downstream, the proper- 
ties of the fluid can be determined if the chamber properties are known. 

MAXIMUM WORK AND EQUILIBRIUM 

Work can be obtained from a system as long as the pressure and 
temperature are not in equilibrium with the surroundings (potential, 
kinetic and other energy forms shall be neglected). Considering the 
system and the surroundings as a unit that cannot receive or lose heat 
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energy, the whole work energy is supplied by the difference between 
' internal energy before (U ') and after (U") the change of state. The internal 

energies of the unit before and after the process are composed of the 
^ internal energies of the system and the surrounding, 

U’ = Ui + Uoi 
U" = U, + Uo2 

where the subscript zero (o) indicates the surrounding. 

Considering that the pressure and temperature of the surrounding 
do not vary, then the surroundings change from the initial to the final 
state by a constant pressure and constant temperature process. At the 
same time, however , the volume of the system changes, and the work done on 
the surroundings at constant pressure amounts to 

Po(Vi - V2) 

Furthermore heat energy (Qq = TqASq) is transferred to the surrounding 
without changing its temperature Tq. The change of internal energy in 
the surrounding amounts to 

U02 - Uoi = ToASo + Po(V, - Vi) 

where ASq = Sq^ - Sq^ indicates the entropy change of the surrounding. 

' The total entropy change equals 

^0 + - Si > 0 
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The work that can be obtained from this unit is 

W < Ui + Uoi - (U2 + U02) 

W < Ui - U2 + (Uoi - U02) 

W < Ui - U2 - TqASo + Po(Vi - V2 ) 

If the process is completely reversible, the maximum work will be 
obtained 

Wmax = Ui - U2 - To(Si - S^j + Po(Vi - V2) ( 13 ) 

Using the following definitions for enthalpy 

H = U + PV ( 14 ) 

and the Gibbs function 

G = H - TS ( 13 ) 

the work equation reduces to = Hj - H2 - To(Si - S2) 

Wmax = Gi-G 2 ( 16 ) 

This equation states an important fact: During such a process, only the 

change of (G) can be converted into work energy while the amount 
(S^ - Sj) is carried away as heat. 

When the state of a system approaches equilibrium, the maximum 
work gradually diminishes and is finally zero. 

Wtnax = 0 = AG (at equilibrium) 
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For a closed system undergoing a chemical reaction, an equilibrium state 
is achieved when the sumiTiation of the Gibb's functions for all present 
compounds equals zero. 

EQUILIBRIUM CONSTANT 

The Gibbs function, or free energy function, is the basis for the 
equilibrium therory. The differentiated Gibbs function yields 

dG = dH - SdT T dS (17) 


dG = dU + P dV + V dP - S dT - T dS 

In a reversible process, with only expansion work involved, the first 
law of thermodynamics yields 

dQ = TdS = dU + PdV 

Substituting dU in the prior equation gives 


dG = V dP - S dT 

For an isothermal process -dT = 0 and the equation reduces to 

dG = VdP 

With the ideal gas law 

PV = n RT 

we find 

dG = dP 


(18) 


The integration between two states results in 

G, - Gi = nRTln^ 

i 1 


or on a mole basis 


G 2 - Gi = RT ln|^ 


(19) 
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This equation can be interpreted as: The maximum work that can be 

obtained in a reversible constant temperature process, going from one 
state with a pressure P ^ to a second state with pressure Pj , is equal to 
the Gibbs function difference. 

Assuming that the Gibbs function (G^) at one atmosphere and all 
temperatures are known, the Gibbs function at any other pressure is 

G = G° + RTlnP (20) 

Formerly it was believed that the reactants in a thermochemical reaction 
were completely consumed to form the products. But developments in 
chemistry have shown that the reactants never disappear. After initiation 
of the reaction, the concentration of the reactants decreases while the 
concentration of the formed products increases correspondingly until 
equilibrium is reached. Guldberg and Waage stated a relationship between 
the concentrations and reaction rates: The rate of reaction is 

proportional to the active concentration of the constituents. 

Consider the following chemical reaction 

aA + + dD 

where A and B are the reactants, C and D the products, and a, b, c, d 
the corresponding mole numbers. After starting the reaction, more and 
more products will be formed from the continuously decreasing reactants. 
During this process the Gibbs function of each species will change 

G = G° + RT In P 
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until a condition is reached in which the same amount of products are 
formed from the reactants as reactants are formed from the products. 

This is the equilibrium condition. It was established earlier that no 
work can be produced from a system when it is in equilibrium with it’s 
surroundings. In the above reaction, the surrounding of the reactants 
is all the products. Therefore, the difference between the Gibbs function 
of the products and the reactants must be zero at this condition, 

AG - 0 = c Gq + d Gp) - (a Gy^ + b Gg) 

Applying equation 20 yields 

0 = cGc° + dGo° - (aGA° + hGe^) + RT (InPc'' + In 
Expressing 


AG° = c Gq° + d Gd° - a G^° - b Gg° 


results in 


( 21 ) 


0 = AG° + RT In 


-p a -p b 


In this equation the following expression is called the equilibrium 
constant, 

K 


( 22 ) 


Since the Gibbs function G is a function of temperature only, Kp 
is independent of pressure and dependent only on temperature. However, 
this does not necessarily mean that the concentrations of the reacting 
components remain the same when the pressure is increased. Since the 
chosen standard condition for the pressure was one atmosphere, the 
equilibrium constant Kp is dependent on the units for pressure. The 
subscript P has been chosen to indicate this. 
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Pressure has an influence on the chemical equilibrium only when a 
reaction occurs with a change in mole numbers or when the gases at 
high pressures cannot be considered ideal. Le Chatelier and Braun dis- 
covered the principle: If a factor influencing the equilibrium changes, 

an effect is produced that tends to reduce the action caused by the 
change. 

CONDENSED PHASES 

A pure substance is defined as having an invariable chemical 
composition for all states of aggregation: solid, liquid and vapor. 

Single-phase or multiple-phase conditions are considered to be in an 
equilibrium state. A pyre substance will have, in general,, the follow- 
ing pressure- temperature and temperature- entropy diagrams (Ref. FIG. 1 

and 2) . 



FIGURE 1. GENERAL PRESSURE- 
T04PERATURE DIAGRAM FOR 
DIFFERENT STATES OF AGGREGATION 



FIGURE 2. GENERAL TEMPERATURE' 
ENTROPY DIAGRAM FOR DIFFERENT 
STATES OF AGGREGATION 
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During a melting (fusion), sublimation and vaporization process, 
the corresponding two-phase conditions are in equilibrium: solid- liquid , 

solid-vapor, liquid-vapor. All three states of aggregation are possible 
only at the triple point. 

In this report condensed and gaseous compounds are treated in- 
dependently, and are expressed by the application of separate equilib- 
rium equations. The vapor pressure is associated with the gaseous 
compound only and is zero for the condensed phase. Therefore, the 
vapor-condensed phase equilibrium is included in this procedure. Solid- 
liquid equilibrium states are treated in Reference 2. 

A reversible isothermal constant pressure process converting 
liquid into vapor can be performed at the equilibrium pressure. The 
entropy change during the vaporization is 



* The change in the Gibbs function is zero because 

^^vap " ^^vap ’ T "^vap = 0 

This can be generalized for any kind of transition between two phases 
of a pure substance, such as solid-vapor or solid- liquid. A pure 
substance is assumed to keep the same thermodynamic properties despite 
the presence of the same matter in another state of aggregation or 
another substance . 

As shown in equation 18, the Gibbs function for a reversible 
isothermal process is 

dG = V dP 

\ 

k\ 
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Performing a vaporization process of a substance, the change in the 

P# 

Gibbs Function is i 

^^vap dP = J (^vap " V^)dP 

^vap 

The volume of a liquid can usually be neglected when compared with the 
volume of the vapor. If, furthermore, the vapor behaves like an ideal 
gas, the equation reduces to 

= J V^apdP = J ^ dP = n RT In 


For a unit mole, the equation reduces to 


AG = RT In 


vap 


In equilibrium, the Gibbs function 
possible when the pressure of the Rquid 


P, = P 


vap 


= P 


equ 


equals zero, 
is the same as 


(23) 

This is only 

the vapor pressure 


Neglecting any effect of the applied pressure on the vapor, which is 
compatible with the assumption that i-S negligible, and treating the 
vapor like a pure substance, the vapor pressure remains a function of 
temperature only. The Gibbs function changes to 


AG^-ap ““ “ ^ equ 


(24) 


Comparison with equation 23 shows that the equilibrium constant is 
equal to the vapor pressure. The pressure of the condensed phase has no 
effect on the equilibrium constant. 
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Example: 


H2 0(i) = 


= - RTlnKp = -RTlnP^ 


AG _ AH - TAS _ AH ^ AS 
RT RT RT R 




^vap Syap 


ni Si 


AH = 10520 cal mole, is the latent heat 


AS = 28.393 cal mole K 


R = 1.987 


mole K 


10520 ^ 28.393 

■ 1.987 X 298 1.987 

Kp = 0.03125 atm 

This result is equal to the vapor pressure at T = 298^K. 

SPECIFIC HEAT AT CONSTANT PRESSURE 
The specific heat at constant pressure is defined as the amount 
of heat added to a specified amount of matter to increase the tempera- 
ture by one degree while the pressure is kept constant. 

c 

■ AT 

Expressing AQ by the first law of thermodyamics 


AQ = AU + PAV 


yields 


AU + PAV 


Applying the definition of enthalpy 


AH = AU + PAV 


we obtain 


Cr. = 
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SPECIFIC GAS CONSTANT 


The equation of state for an ideal gas is 

Pv = RT 

At constant temperature and pressure, one mole of two different gases 
may be described as 

Pvx Mx = Mx Rx T 

Pv^ M2 = M2 Rz T 

Avogadro*s law states that the mole volumina (v-M) of different gases 
are the same when the temperatures and pressures are equal. 

VxMx = V2M2 


or 

Ml Rx = M2R2 = MR = R 


(27) 


R is called the universal gas constant. 

For a mixture of ideal gases enclosed in a volume (V) j the total 
pressure is equal to the summation of the partial pressures of all 
constituents when the temperature is homogeneous. 


’tot = ? Pi 
1 


(28) 


Replacing the pressures by the corresponding term of the general gas 
equation, the specific gas constant of the mixture can be calculated. 
Mtot Rtot ^ _ y Rj 

17- “ 2^ V 


SMiRi 

^tot = 


( 29 ) 


tot 
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Since the mass is equal to the product of the mole number and the 
molecular weight, 

m = nM (30) 


the specific gas constant can also be determined without knowing the 


specific gas constant of each constituent. 

Rtot = 

^ mtot 


R = MiRi (31) 

2 Hi _ 2 

mtot S nj Mi 

If a condensed phase is existent, its share is only considered in the 
total mass, not in the total mole number, because we neglected the 
volume as previously indicated. 

DALTON’S LAW 


Dalton’s law states that the total pressure of a gas mixture is 
composed of the partial pressure of the individual components. 


Ptot=fPi 

Hereby, a gaseous constituent fills the whole volume and behaves as if 
the other gases were not present. When the temperature of all the 
gases is the same, 

Pi V = m'l Ri T 


P^V = m 2 R 2 T 


Ri 

P2 mj Rz 

Since the mass equals the mole number times the molecular weight, 
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we obtain 


m = nM 

Pj _ TLi Mx Ri 

P 2 112 M2 R2 

Using the relationship between the specific gas constant and the 
universal gas constant 

R = MR 

Zi = ill 

P2 ~ 1^2 
or 

P _ n 
^tot ^tot 

The pressures can also be expressed by the general gas equation 

P - nRT 
V 

Defining the concentration to be 

n 

^ V 

we obtain 


P = CRT ( 32 ) 

This relationship can be used to express the equilibrium constants by 
concentrations or mole numbers. 
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HEAT OF FORMATION 


When a compound is formed from reactants, the heat liberated or 
induced is called heat of formation and represents an enthalpy change 
during the reaction. 

Mostly, thermodynamic processes occur in fixed nonreactive gas 
mixtures. In this case, use can be made of the existing thermodynamic 
property tables for the given substances, since only changes in state 
of these substances must be considered. The enthalpy tables are based 
on a reference level that is arbitrary since only enthalpy changes are 
signi f icant . 

When a chemical reaction takes place, this concept is not adequate. 
Compounds may change or new ones may be formed from the reactants. When 
an energy balance is made between the reactants and reaction products, a 
common reference level has to be established. Since reactants and 
reaction products are composed of the same elements, the enthalpy 
reference level for the elements is arbitrarily chosen to be zero at 
standard conditions (1 atm pressure and 298.15^K). The enthalpy of 
formation determines the enthalpy level of a specific compound at 
standard conditions with respect to the chosen reference level. If the 
heat of formations are summed for the reactants and products separately, 
their difference indicates how much energy is liberated or has to be 
induced during a chemical reaction. 

In case of a combustion reaction, the energy difference is 
available to heat the combustion products from standard conditions to 
a higher temperature. 
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STOCHIOMETRY 


By definition a stochiometric mixture ratio corresponds to complete 
combustion, the formation of saturated molecules in the combustion 
product. A method to solve this problem is based on the principles of 
oxidation and reduction. 

Definitions: 

Oxidation occurs when electrons are lost. 

Reduction occurs when electrons are gained. 

The loss of electrons is represented by a postive valence; the 
gain of electrons is represented by a negative valence. The valence, 
therefore, determines the degree of combining power of an element or a 
radical. 

General rules (Ref. 7 and 8) 

1. Free elements have no valences. 

2. Hydrogen in combination always has a positive valence. 

3. Oxygen in combination has two negative valences. 

4. The sum of the valences in any compound equals zero. 

In an oxidation-reduction reaction, the total increase of positive 
valences equals the total increase in negative valences. 

Example: Formation of water 

Hz + i 02 - H2 0 

including the valences 

+ i 02° 0-2 
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In this reaction, all compounds have a total valence of zero; however, 
in H^O, the valence of hydrogen is (*f2) and the valence for oxygen is 
(-2). Hydrogen, therefore, has encountered a loss of electrons (oxi- 
dation), and oxygen, an increase of electrons (reduction). 

The atoms of the various elements differ in their combining capaci- 
ties and, therefore, have different valences. To determine the stochio- 
metric mixture ratio for a combustion process from the reactants, it is 
assumed that in all formed compounds, the elements possess a valence 
that is the most common. 


element 

C 

F 

H 

HE 

N 

0 


valence 
+ 4 

-1 

+1 

0 

0 

-2 


If a reaction takes place between oxygen and hydrogen and the 
products formed possess the most common valences of these two elements, 
we can transfer the valences to the reactants, because all products are 
composed of the reactants. 

V V 

TT F 0 

Products 

V = valence 

^ = atomic number of an element in the fuel 
^ = atomic number of an element in the oxidizer 
M= molecular weight 
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In one gram of hydrogen there are 


valences • 


Y 

In one gram of oxygen there are ^ ^ valences. 

0 Mq 

Since the valences of oxidizer and fuel in a stochiometric com 
bustion must be equal, we can form a ratio of the two values, including 
a proportionality factor, so that this ratio equals one. 

- 

-(vp f p Me) * ® 


The minus sign makes the ratio independent of the sign. 

Solving for the stochimetric mixture ratio (4>s), we obtain 

vp Me 
‘♦’s " - ve^e^^F 


When each reactant is composed of more 
must read: 

(s viei)^Me 


elements , 


as in CjHj, the equation 

( 33 ) 


If one of the propellants is a mixture of two or more compounds, as 
Oj ^ ^2’ summation of the parenthetical terms must consider the 
percentage quantities of each compound. 
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COMBUSTION AND EXPANSION PROCESS IN A ROCKET ENGINE 


Most liquid propellant rocket engines use two propellants, a fuel 
and an oxidizer. Each propellant, is forced into the combustion 
chamber by a pump or by a high-pressure syscem. Nozzle walls consist of 
many small tubes brazed together, and are regeneratively cooled by one 
of the propellants, recovering most of the transferred heat energy and 
feeding it back into the combustion chamber. Another important device, 
the injector plate, consists of a series of small orifices that mix the 
fuel and oxidizer in the forward end of the combustion chamber. The 
quality of mixing influences the space occupied by the combustion re- 
action and, therefore, is important in determining combustion chamber 
length. Combustion stability is also affected by the quality of the 
mixture. 

Cryogenic propellants usually are injected at the boiling point, 
unless one of the propellants is a mixture (for example with the 
oxidizer liquid fluorine - liquid oxygen, the injection temperature is 
the fluorine boiling point and the liquid oxygen is undercooled). 

Storable propellants are normally injected at ambient temperatures. 

The isenthalpic reaction of the fuel and oxidizer liberates heat 
and produces combustion products that dissociate and recombine. When 
the mean composition of these products stabilizes (that is, the rates of 
recombination and dissociation are equal), a state of equilibrium is reached. 
The potential energy of the combustion products is isentropical ly con- 
verted into kinetic energy when expanded through the nozzle. A choked 
condition develops at the throat when the following equation is satisfied: 
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( 34 ) 


1 



At this point, the velocity of the gas particles equals the local speed 
of sound. Since this velocity can only exist in the throat area, the 
niass flowrate can be calculated. 

The combustion of the reactants does not occur instantaneously in 
a specified area in the* combustion chamber but rather in some portion of 
the chamber. During this period, the quantity of gases formed progresses; 
this naturally requires an increase in velocity. The velocity increase 
calls for a drop in static pressure according to Bernoulli's equation 
for compressible flow (Ref, 9). 

2 Poo/" 

Pst - Poo = Woo' — + 

(The equation is only valid in the range of Ma - 0 to Ma = 0,9) 

The continuing combustion, therefore, releases energy at constantly 
falling pressure levels. The total energy released during the reaction 
is a fixed value, but the availability of this energy varies. The com- 
bustion occurring in a finite length causes a pressure loss and an in- 
crease in entropy. 

With rising temperatures, the reaction rates of dissociation and 
recombination increase. Therefore, at high temperatures a change of 
state of a system reaches equilibrium quickly; however, at low tempera- 
tures more time is needed to reach equilibrium. 
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During the expansion process, the combustion products change 
state continuously. When the reaction rates are so fast that a chemical 
and thermal equilibrium is maintained at any point, the condition is that 
of shifting equilibrium. If the reaction rates are slow, freezing 
almost any dissociation or recombination during the expansion process, 
the condition is frozen equilibrium. 

THEORETICAL COMBUSTION AND EXPANSION MODEL 

The following assumptions are made: 

1. Propellants are injected at boiling conditions or at ambient 
temperature. Injection velocity is negligible. 

2. Combustion occurs instantly in an isenthalpic process at 
combustion chamber end stagnation pressure and a predetermined mixture 
ratio. 

3. The expansion process is isentropic and follows either frozen 
or shifting equilibrium conditions by minimizing the Gibbs function. 

4. All gases behave according to the ideal gas law. 

5. Condensed phases and the corresponding gas phase possess 
the vapor pressure. 

6. The occupied volume of a condensed phase is negligible. 

7. No heat transfer occurs through the nozzle wall, 

8. Thermodynamic and performance parameters are uniform in a 
cross-sectional area of the nozzle. 

9. The velocity of the gas and condensed products are the same 
and are parallel to the nozzle axis (one dimensional). 
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W II 


10. No boundary layer influence is considered. 

11. The performance parameters are a function of expansion ratio, 
but are independent of nozzle wall contour. 

12. There is no jet separation from the wall towards the nozzle 

exi t . 

DEVELOPMENT OF THE EQUATIONS USED IN THE 
COMBUSTION AND EXPANSION CALCUT.ATIOM 

In a combustion process at constant pressure, compounds formed 

from the reactants can be represented by the following equation 

^bF' ^cF • • • ) + ^6 . . .) - 

100 110 115 

JT'i (Zai Ybi ^ci • • • ) + I nj ( Z^i Ybi Xd . . . ) + Z ni(ZaiYbiXci 
1 i = 101 i = 1 1 1 

where 
100 

^ ni(ZaiYbiXci ...) represent the gaseous compounds 
i= 1 

1 10 

(^ai ^bi ^ci •••) rep sent the gaseous elements 

i - 1 0 1 
115 

Z (^ai '^bi ^ci . . . )i'epresent the condensed products 

i = 1 1 1 

All equations involve the mole numbers and the partial pressures. 
To simplify the problem, the chemical reaction equation is multiplied 


(36) 

...) 
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by a factor (A), which makes both the mole number and the partial 
pressure of each compound numerically equal. 


A[Np (Zap XcF ...) + Nq ( Za q q 


...)] = 


( 37 ) 


100 110 115 

( ^ai ^bi ^ci •••)+ ^ (^ai^bi^ci ••») 

1 i = 101 .1=111 


A = Nj 

The equilibrium state can be determined by the following equations, 
1. Chemical equilibrium equations 

Consider that each reaction product is formed from gaseous 

atoms . 

Z + bi Y + ci X ~ Zr^^i Y^i X^i 


The equilibrium condition for these reactions follows from equation 21 


AGi = AGi° + RT (In Pi - ailnPz - bi In Py - ci In Px - • • • ) (38) 


AG| must equal zero at equilibrium. For any other condition close to 
equilibrium, will be a small number 6^ so that 

5. = ( ^ ) + In Pi - ai In Pz - bi In Py - In Px - ... (39) 

VRT/i 

For condensed products, the In terms are zero because they have 
no influence on the equilibrium constants. 


A3 




6i = :r— - ai In Pz - bi In Py - Ci In Px - • • 

\RTji 


(40) 
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2 . 


Mass balance equations 

The total mass does not change in a reaction; also, the 


masses of the individual elements stay the same. 

1 

^tot = A ^ 
i = 1 


These masses 
reactants. 


1 

^ot = T Z biNi 
^ i=l 


1 

*^tot ” A ^ 


i = 1 


mu 


St be equal to the masses of the elements in the 


(41) 


3. Pressure equation 

In a constant pressure reaction, the total pressure remains 
the same. Due to Dalton's law, the summation of the partial pressures 
of all compounds in the reaction product must be equal to this pressure 

no 

P = Z Pi (42) 

i= 1 

4. Constant enthalpy equation. 

As was mentioned in the section ’’Equilibrium Constant", the 
enthalpy of the reactants must equal the enthalpy of the combustion 
products. The enthalpy of reactants is equal to the heat of formation 
plus the heat energy required to go from the reference level to the 
injection conditions. 
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Cp°dT)eNe 


(43) 


HO = S 
F 


± -L 

(H£ + f Gp°dT)pNp + 2 (Hf + f 


ref 


ref 


The enthalpy of the combustion products is equal to the sum of the 
individual enthalpies of the formed compounds. 


H=is(HT»)iNi 


(44) 


Note: The heat of formation is already included in the terms. 

5 . Constant entropy equation. 

The entropy of the combustion products is given by the 
following equation as developed in the section "Entropy". 

For gaseous compounds 

St =;^ r[(ST'’)i - RlnPilNi (^5) 

For condense,d products 

St 2(ST°)iN£ (46) 

To determine the equilibrium condition for the combustion process 
at constant pressure, the equilibrium, mass balance, pressure, and 
enthalpy equations must be solved simultaneously. During the expansion 
process, the same equations are considered to determine equilibrium 
conditions except the enthalpy equation is replaced by the entropy equation. 
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Both sets of simultaneous equations are non-linear. The Newton-Raphson 
method is applied to obtain a solution using a finite difference 
approximation to the total differential. 


NEWTON-RAPHSON METHOD 


Consider two nonlinear simultaneous equations in x and y in the 
following form 

Fj = Fi (x, y) = 0 
Fj = F 2 (x, y) = 0 

When X and y are the solutions there exists a difference for any 
other estimate of x^ and y ^ so that 

Fi (x, y) - Fi (xi, yi) = AFi if 0 


F2 (x, y) - F2 (xj, yi) = AF2 ^ 0 

The definition of the difference is arbitrary (as far as signs are 
concerned), but when it is once established, the proper correction 
equations must be set up. 

The total differentials of the given functions are 


9Fi , or I . 


9F, 


dF, = 


9F, 


dx + 


dF, 


dy 


9x ■■ 9y 

Going to finite differences, the equations change to 


AFi = 
AF2 = 


(s)“ • (t)' 


(47) 
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These equations are linear in Ax and Ay when numerically solved at the 
point x^and y^^. Because of the transition from the total differential 
to the finite difference form, the equations are only approximate. 
Therefore, the correction variables are 

Xi+ 1 = + Ax (48) 

yi + 1 = Yi + 

+1 + 1 ^PP^o^^tiing X and y. 

Using the following mathematical concept 

d In X _ 1 
dx X 

or in finite difference form 
A In X 1 

. 'V 

Ax X 

Ax ~ X Alnx (49) 

the total differentials can be rewritten. 

According to these equations, the set of nonlinear simultaneous 
equations is transformed into a set of linear simultaneous equations. 
The corresponding correction equations are 

Inxj^ 2 = Inx^ + Alnx 01) 

In Yi ^ 1 = In yj + Ain y ' 
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This procedure is used to solve for the equilibrium condition in the 
combustion and expansion process. 


RESULT OF EQUATIONS LINEARIZED BY THE NEWTON-RAPHSON METHOD 
The final equations, as used in the equilibrium calculation 
solving for N^ = , T, A, are presented below. A matrix can be 

constructed as shown on page 34 (FIG# 3)# The equations are derived in 
Appendix B . 

Equilibrium equation for gaseous compounds 

AF = Ain Pi - ai AlnPz - biAlnPy - ci AlnPx - . . . - [— — ^InT (52) 

\ Ji 

Equilibrium equation for condensed phases 

AlnT (53) 

Mass balance equations 

115 115 

AFa = - -7- T aiNiAlnA + Y T aiNiAlnNi (54) 

1=1 1=1 

Total pressure equation 


f^T 

AF = - ai Ain Pz - bi Ain Py - ci Ain Px - . . . - I -= — 

y RT 


no 

^P= S Pi Ain Pi 
i = 1 


(55) 
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Enthalpy equation 

115 _ 115 _ 

^H=-X 1 (HT°)iNiAlnA + - Y, (HT°)i Ni Ain Ni 


A '“-L “ ■ A . , 

1=1 1=1 

115 _ 

+ T Y (Cp°)iNiAlnT 

i= 1 


Entropy equation 


a) for gaseous products 
1 110 _ 

^S=-X T. [(ST°)i - RlnPilNiAlnA 

■ 110 _ 110 _ 

Y (ST°)iNi AlnNi - Y RNilnP^AlnN^ 

.1=1 i= 1 

no _ 

- Y ^in Pi 

i= 1 

b) for condensed phases 




1 - 

+ 4- Z <Cp")iNiAlnT 


1 = 1 


( 56 ) 


(57) 


115 _ 115 

^S = -T ^ (ST°)iNiAlnA + - Y (St“)i Ni Ain 
1=111 ^ . 


1=111 


115 


+ T Z (Cp°)iNiAlnT 


1=111 


(58) 


Combining both equations yields 

115 110_ 110_ 

Y (S-p°)iN. AlnN. - Y RN.lnP.AlnN.- Y RN^AlnP- 

1= 1 i= 1 


AFo = 7 - 

S A 


+ T Z (CpO)iNiAlnT - 


115 


AFs =- 


1 

A 


1=1 

115 no _ 

Z (St°)i- Z 
1=1 1=1 


C'l: 

i+ DjNi' 


1 = 1 

_ 110 _ 
(ST°)lNi - Y RNilnPi 


1 = 1 


Ain A 


115 


iAlnPi + ;^ Z (Cp°)iNiAlnT 


1 = 1 


A 


115 110 _ 

Z (Sj )lNi - Z RNilnPijAlnA 


1 = 1 


1 = 1 
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FIGURE 3. ORIGINAL COEFFICIENT MATRIX FOR THE SET OF SIMULTANEOUS 
EQUATIONS SOLVING FOR Ni, A, AND T 




FIGURE 4. REDUCED AUGMENTED MATRIX USED TO SOLVE 
FOR Ni, A, AND T SIMULTANEOUSLY 














DERIVATION OF THE MOLAR SPECIFIC HEAT AT CONSTANT PRESSURE 


By definition: 

AdH = NCpdT 

(59) 

A dH 
N dT 


The total enthalpy of the product is equal to the summation of the 
proper enthalpy portions of the compounds. 

115 _ 

H=t T (HT°)iNi (60) 

^ i=l 

The enthalpy is a function of: 

H = F(A, N, Ht°) 


At first, the enthalpy equation is differentiated with respect to T in 
order to obtain an expression for the Cp equation. 

1. In the case of shifting equilibrium, we obtain; 


8H - 9H — 0. 9H _ 

dH = dNi + 0^ d (H-j- )i + ~ dA 


9N; 




The individual partial derivatives are 


a A 115 \ , 115 

i = 3 (Sa(a Ni 


8H 




9H _ 8 
8A ■ 8A 


A 115 _ \ 1 115 _ 

^ Z (HT°)iNi) = 1 (HT°)iNi 


Substituting the proper terms and dividing by dT yields. 

115 _ 115 _ 1 115 _ 

I (HT“)idNi+;j I Nid(Hx°)i-;^ I (HT'UNjdA 

i=l i=l i=l 


115 


115 


. dT " A )i dT ^ A dT “ A^ )iNi 

1=1 1=1 1=1 

According to the mathematical concept 

d Inx _ 1 
dx X 


and with 


we obtain: 
115 


d(HT°)i _ _i.f fc 


dT 


^t) i - (Cp )i 


dH 


dT A . 


115 


115 


= T (Cp'UNi-^ 2 (hANi 


i= 1 


i= 1 


AT . 


i= 1 


din A 
dlnT 


dH 1 
dT " AT 


115 115 _ 115 _ 

I T Z (Cp'UNi - Z (HT°)iNi 


d In A 
dlnT 


. , dinT . . , 

1=1 1=1 1=1 

The molar specific heat at constant pressure is now 

= Til— (f ) = -n^[ .'i:' <«T”)i Ni ^ ^ T 'z (Cp-)! Nj 

V V L‘ = l ‘=1 


Z N, 

i= 1 


T I Nj 
i = 1 


i = 1 

Since Ni = Pi, the term ^ ^1 can be replaced by the 

d In T d In T 

corresponding expression obtained from the equilibrium equation 52. 


37 



( + inPi - (ailnPz + bilnPy + cilnPx +...) = 0 

V RT ji 

dlnPi__ dlnPz ■ , dlnPy . dlnPx x 

__ = ai +Ci j 


dlnT dlnT dlnT ’ dlnT ’ \ 

The molar specific heat at constant pressure finally results in 

Cp 


1 _ « r dlnP -7 , dlnPv , dlnPx . 

= -Yi| \,l (HtUNi U i-j^ + bi-3^ + ci-jj^+... 

V TVT. ^ ^ 


T X Ni 
i = l 


(62) 


2. For frozeti equilibrium and A, are no longer functions of 
temperature. Therefore, equation 62 reduces to 


Cp = 


1 r — 

ns— T I (Cp'UNi 

V N. L 4=1 


T X Ni 
i = 1 


115 


I (Cp°)iNi 


(63) 


C - 

^P~ 115 


Z Ni 
i= 1 


DERIVATION OF THE MOLAR SPECIFIC HEAT AT CONSTANT VOLUME 
Applying equations 156 and 157 to equation 118 in Appendix B yields the 
following equation for Cy 

(li)§ (6^) 

7W) 

V8P/t 


C., = C, 
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Take the difference of the specific heats from equations 156 and 157 


Cp 



(65) 


The derivatives are calculated from the property equations: 
Combining both equations 140 and 141 yields 


dT 


(9t)p ■ (st)^ = (avX ■ (ap)T 


Comparing this with equation 65 shows 


dT 

T 


Applying equations 144, 146, 147, 149, 152, 153, results in 


dT = 


Cp - C, 


(It) 


dv + 77- 7 “ 


i \ 

\8T/] 


dP 


( 66 ) 


Comparison of the coefficients with equations 142 yields 


/8t\ t / ap \ 

y8v /p Cp - Cy\8T/^ 

/8T\ T / Bv \ 

\BP)y ^ Cp - \BT)p 


Both equations have the same result 



from equation 118 gives 


Substituting 



V 


:p - c^= 


(S)‘p( 


8P \ 
8v )r^ 


( 67 ) 


{li)r 


m. 


In the following section expressions are developed for 

Starting with the equation of state for an ideal gas 

MPv = RT 


and 


taking the logarithm, forming the differential, and dividing by dlnX 


yields 


Since 


d In M d In P din v _ dlnR ^ dlnT 
d In T d In T d In T " d In T d In T 


dlnR 


dlnT 


dlnT 

dlnT ■ 


It follows that for a constant pressure process 


/ d In V \ , ( dlnM\ 

\^dlnT/p " ■ ^dlnT/p 


/ dv \ . dlnM~ | (68) 

\dT/p " T L ■ dlnTjp 
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Beginning again with the general gas equation, taking the logarithm, 
forming the differential, and dividing by din P yields 


dlnM ^ dlnP ^ dlnv _ d In R dlnT 


with 


d In P d In P d In P d In P d In P 
dlnR 


dlnP 

dlnP 


= 0 


= 1 


dlnP 

For a constant temperature process 


or 


/ din V \ _ , / d lnlvl\ 

ydlnPy-j- ydlnPy-p 

/ dv \ _ V d In m1 

\dP)rj, " “ pL^ dlnPj-p 


(69) 


Solving equation 67 for Cy and substituting equations 68 and 69 for the 
corresponding terms results in 


Cy = Cp - 



dlnM 

2 

vP 

^ ■ dlnT 

P 

T 

‘ dlnM' 

^ d In P 



T 


(70) 


Changing over to molar specific heats according to 

Cy = MCy 

Cp = MCp 

and using equation 70, we finally obtain 


(71) 


Cy = Cp - R 


1 


d InM 


dlnT 


1 + 


d In M 
dlnP 


2 

P 


(72) 
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RESULTS OF EQUATIONS LINEARIZED BY THE NEWTON-RAPHSON METHOD TO SOLVE 
FOR THE PARTIAL DERIVATIVES AT CONSTANT PRESSURE AND TEMPERATURE 

As shown in Appendix A, the equation of state can be expressed by 


F = F(P, V, T) = 0 

-mMMi 

Since any other property can be expressed by the three properties, P, v, 
T, any other partial derivative can be expressed by the three partial 


derivatives. 

In this report, the following first partial derivatives are used 



/ d lnlvl\ 

ydlnT /p 

/ dlnM \ 

y d In P / ^ 


This leads to the determination of the specific heat at ’constant pressure 
Cp and constant volume Cy These terms can be expressed by the quantities 



To solve for the partial derivatives, a set of equations must be 
solved simultaneously. The soultions of the equations, which are 
derived in Appendix C, are given below: 

1. The derivatives at constant pressure can be obtained from the 
following equations: 


♦ 


♦ 
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a. Equilibrium equation for gaseous compounds 

dlnPi dlnPz dlnPy _ _ 

b. Equilibrium equations for condensed phases 

^ /aHt°\ dlnPz V. llilPX . - ... 

= V'rT’A " ^i'dhTr ■ dlnT dlnT 


c. Mass balance equation 

dlnNi 


“5 dlnA " = 


^ ainA , ^ i^T. uui IN] 

0=- Z dlnT dlnT 

i= 1 1=1 


( 73 ) 


(74) 


(75) 


d. Total pressure equation 

0 = p. 

^ L ^1 fllTI T 


i= 1 


(76) 


This simultaneous set of equations can be simplified when the equilibrium 

equations for gaseous compounds is solved for *1 1*^ l*i and the result is 

d In T 

substituted in the mass balance equation for 

d In T 


This yields 

115 


- dlnA , 

0 = - Y. ai Ni -T7:-7r + 


i = 1 


dlnT 


110 r 


dlnPz j ^ dlnPy 
dlnT °i“TTr^ 


dlnT 


V RT /ij >1 > Id 


'1 dlnT 


InNi 

dlnT 


(77) 


In writing a coefficient matrix for these equations, using equations 
77, 74, 76, the system is incorporated in the reduced augmented matrix 
solving for the mole numbers. 

2, The derivatives at constant temperature can be obtained form 
the following equations: 
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dh^ 

din A 


b. 


Equilibrium equations for gaseous compounds 


dlnPz , dlnPy _ dlnPx 
dlnA ■ dlnA ' dlnA 

Equilibrium equations for condensed phases 


(78) 


n = dlnPz , dlnPy ^InPx 
^ dlnA " i dlnA ^ dlnA 

c. Mass balance equations 

U5 „ 

0 = -aoA+ Y, aiNi^^i2Lfi 


(79) 


i = 1 


d In A 


(80) 


This simultaneous set of equations can be simplified by replacing 
d In P. 


the term 

d In A 

equation. The result is 

110 

0 = - aoA + 5) aiNi 

i = 1 


Jl with the proper expression from the equilibrium 


a- + c- + 

^dlnA ^dlnA ’■dlnA 


(81) 


115 „ 

I ^ -NT ^ Pi 

+ Y ^ 

i = 1 


^ dlnA 


The coefficient matrix for this simultaneous set, using equations 79, 
81, is also incorporated in the reduced matrix. 

The relationship between ( ^ ) and ( ] is derived in 

^ \d In T/p \d In P/.p 

Appendix D, 

ISENTROPIC EXPONENT 
Speed of sound is by definition 




M 


(82) 


44 



Using the adiabatic relationship for a perfect gas, 

= P — ^ = const (83) 

taking the logarithm 

lnP-Ylnp= const 

and forming the differential yields 
dlnP - ydlhp = 0 


dlriP _ pdP 
■ dlnp ~ Pdp 


With the relationship 


we obtain 



V 


lnp= In 1 


In V 


In 1 = 0 


(84) 


d In p = - d In V 
dp _ dv 

p " " V 

such that 


vdP 
■ P dv 


Setting both terms for Y equal to each other and solving for 


d P 
d P 


yields 


dP _ PvdP 
dp “ " pPdv 


Substituting this in equation 82 gives 


Wc 


/ Ei(^\ - ITil^S 

V" /s ■ P \dv ) 


RT 


(85) 
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The specific heat ratio is by definition 

Cp 

^ = 7=^ 

L.V 

from equation 64, it follows that 


/8v \ /8P\ 

Cv " 


Applying equation 69 results in 


dlnM l /8P\ 

~ ^ L ^ dlnP Jx\9v yj 

_ Fi ^ d'lnMl 
dlnP 

v_/^\ 

5 ■ P\9v/s 


This expression appears in the equation for the speed of sound. Rep lac- 


te). 


by equation 64 results in 


" p 




Considering furthermore equations 67, 68, 69 gives 


= ■ p 


V 

Vs ” TD 




S "p Cpvf, . dlnMl . Tv^ [, dlnM ' 
_~^P~ ^■^dlnPU dlnT 


Applying the general gas equation and going over to molar specific heats 
finally yields 


r. dlnMl dlnMl^ 



THERMODYNAMIC PROPERTIES 


The thermodynaraec properties per mole of the individual combustion 
products used in the calculation procedure are the specific heat at 
constant pressure Cp, the enthalpy and the entropy S^. These 
properties are correlated by the following definitions: 

Cp° = F(T) 

Ht" = / Cp® dT 

s” dT 

When the Cp values are represented by the following polynomial 

Cp® = a + bT + cT^ + dT^ + eT^ (88) 

in which a, b, c, d, e are constant values, the enthalpy and entropy 
equations yield correspondingly 

Ht® = aT + I + I + I + f T® + const (89) 

St° = alnT + bT + -| + I + -| + const (90) 

The enthalpy function must include the reference enthalpy H°298.15 
298.15 K. In most cases, the properties are tabulated as a function of 
temperature at atmospheric pressure. To minimize the errors that show 
up if each polynomial is developed individually, all three polynomials 


i 
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are developed simultaneously by the least-squares method. This is 
advantageous because all corresponding coefficients have the same value, 
and it saves storage space in a computer. One polynomial does not cover 
the total range from 300 to 5000°K with sufficient accuracy; therefore, 
two ranges were considered, one from 300 to 1000*^K, and another from 
1000 to 5000°K. In this case, however, a constraint must be placed on 
the development of the polynomials so that at 1000*^K the intersection 
point of the two corresponding polynomials have the same value (but not 
the same tangent). 

The equations in the calculation procedures use the properties 


and the polynomials were developed as 
C-o° 

-r=- = A + BT + CT^ + DT^ + ET"^ 


Ht° BT CT^ . DT^ ^ ET^ ^ F 

^ = A + ^ + — + ^ + — +T 


St CT^ DT^ , ET^ , ^ 

-Jr- = A In T + BT + -Y~ + -J- + 

R 


The polynomial coefficients A, B, C, D, E, F, for various compounds and 
elements have been published in Ref. 3. 
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THROAT AREA CONDITION 


Knowledge of flow and thermodynamic properties, of the combustion 
products in the throat area is necessary to determine performance para- 
meters for a rocket engine; for example, characteristic velocity, flow 
rate, expansion ratio, thrust, etc. The approach to compute the throat 
condition is iterative and is based on the criterion that the flow 
velocity is equal to the speed of sound 


w = w. 


Ma = = 1 


(94) 


With known chamber data, a first approximation of pressure in the throat 
region can be made according to the thermodynamic equation derived for 
ideal gases and constant isentropic exponent. 


" ^ch 



(95) 


The corresponding temperature is 

T. = (96) 

The velocity and the speed of sound must be calculated for this condition 
considering a constant entropy expansion process. If the results do not 
satisfy the criterion w = Wg a new approximation of the pressure will be 
performed using the following concept. 

Every thermodynamic property can be determined when two other properties 
are known. 
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Choosing 


h = F(P, s) 

which reduces to 

h = F(P) 

for a constant entropy process. From experience it is known that there 
is a logarithmic relationship between the pressure and the enthalpy. 
Differentiating with respect to In P results in 

or considering finite differences 

During an isentropic process the total energy cannot change. However, 
as long as the proper static pressure in the throat area has not been 
found, the difference between the chamber enthalpy and the total enthalpy 
for the throat region is not zero. 


= hch - ^ 

Applying the mathematical relationship 

AP Pk+i- Pk 
AlnP pj 

with k indicating the previous calculation of the static pressure, we 


obtain 


= 



Pk+1 - Pk 

Pk 


Solving for P ^ yields 


Pk + 1 = Pk 



( 97 ) 
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The total enthalpy in the nozzle is composed of a static and a dynamic 


portion 


^ ~ ^stat l^dyn 

h = hsut t 


(98) 


Multiplying and dividing the dynamic term by Wg^ and the molecular 
weight gives 


h = hgtat + 


Ma^ ws^ M 

2g„M 


Expressing the speed of sound Wg by the following thermodynamic relation- 
ship 


results in 


h = hgtat + 


Ma^ 

2M 


\s RT 


This equation must be differentiated with respect to In P to adapt it to 
equation 97. 


/ 3h \ _/8hstat\ . , Tr rVM) ] 

\ainP/s 81nP /s 2 ^s^\ainP/g 


(99) 


Considering the Mach-number and the specific heat ratio as independent 
of the pressure, the derivative of the dynamic term is derived. From 
the ideal gas law the following expressions are developed 


PvM = RT 
1 

V = — 

P 
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M -.f!©! «,i 

L 9P Js "rL apj ■ rL J"rdL papj 


Applying equation 84 


results in 


/ 9 In p \ _ Jl_ 

\ a In p/s ~ Ys 

M -lI 

I 8P J. RdL 


\M/ ^ 1 Ys - 1 

L 8P Js RpL Ys J 

Next, the static enthalpy derivative is determined. The first law of 

thermodynamics states: 

^stat = u + Pv 


a^stat 

ap 


u 


ahstat 
a InP 


) =- 

Is P 


Substituting the corresponding terms into equation 99 yields 


(alnp), “[ 


Ma 

1 + — (Ys - 1) 


Applying the ideal gas law yields 

PM 

= RT 

P 


/ ah \ RT r 
VainP^g = M L 


1 4. Ma% , , 

1 + — (Ys - 1) 


At the throat area, the Mach-number equals one. Therefore, the last 


equation reduces to 


/ ah \ ^ 

VainP^s ■ 2M 


(Ys - 1) 


(102) 
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With this expression, the following approximations for the static 
pressure at the throat area can be made. 


Pk + 1 = Pk I 1 + 


hch - hgtat + 

1^1 

2M 

+ 

1 


X 


(103) 


EXIT AREA CONDITION 


To determine the condition at a given expansion ratio € , the 
computation is similar to that used for the throat condition. At first, 
a proper static pressure must be estimated for the expansion ratio 
being investigated. An equation from Ref. 4 serves as a basis for con- 
sidering the flow of an ideal gas through an isentropic nozzle. 



(104) 


Since the exit static pressure cannot be determined implicity from this 
equation, the following procedure was applied: 

p 

1. For any combination of the pressure ratios 0.1; 0.01; 

0.001; 0.0001; 0.00001, and specific heat ratios Yg = 1.10, 
1.15, 1.20, 1.25, 1.30, 1.35, 1.40, the expansion ratio was 


calculated. 

P A 

2. Plotting the data versus — — for constant specific heat 

ratios on log - graph paper indicates a fairly linear relation- 
ship. 
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Pch 



FIGURE 5. PRESSURE RATIO VERSUS EXPANSION RATIO 
' FOR CONSTANT SPECIFIC HEAT RATIOS | 

A range of 3 to 100 for the expansion ratio of all specific 

heat ratios’ (Yg) was taken into account. 

Because of the linear relationship of a constant specific 
heat ratio curve, an equation of the following straight line 
equation can be written for each case. 


y^s = ^ 


or if log-scales are used 




+ Inc 


The slope of each constant Y-line can be calculated from the 
end points. 


Knowing the slope and considering one point on the correspond- 
ing Yg-line, the term In c can be computed. 

Plotting the different slopes (m) and constants (In c) as a 

function of the specific heat ratios ( Yg) , a linear relation- 
ship was again found so that a straight line equation could 
be written for either case. 


my - 

’si ^S2 


Inc.^ - In Cm 

(y. - V,.) 

*S1 ^S2 


A general equation can be developed that represents the family 
of curves 


y = m(y )x + b(\g) 


In: 


ch 


'/my - my \ 
[__Isl 

V Vsi - Vs2 / 

{( 


SI ^ „ 

>sl ^S2 
In Cv - In C: 


In 


At 


^S1 " ^S2 




SI + 


or 


^ch 


=: e 


'mv - 
^si 


^S1 ^S2 




- vs,) + 


.] 


■( 


In Cm - In Cm ^ 

»S1 

^si ~ "^sz i 


(\s - 


si 


(105) 


55 




56 









FIGURE 7. ORDINATE SECTION OF CONSTANT 7 LINES 
VERSUS SPECIFIC HEAT RATIO ® 
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9 . 


P 

Pch 


Using real values from Figures 5, 6, 7, the equation reads 

- [{0.867>/g + 0.27165) ln^+ (0.699228 ys + 0.442199)] qq6) 


10, From this equation an approximation of the static pressure 
(P) at a desired expansion ratio can be made. 

With the approximation of the static pressure in the exit 
area, the desired expansion ratio ( e ) will not be obtained 
the first time. Consecutive approximations of the static 
pressure can be calculated from the generally developed 
equation 97 by introduction of an expansion ratio relationship. 

The total enthalpy does not change during the expansion process, 
and equation 98 can be expressed by the expansion ratio (e ). 

Since the mass flow rate through the throat section is equal to 
the mass flow rate through every other cross sectional area, the follow- 
ing equations result : 


At Pt^t = A p w 


e = ~ = ^tPt 
At w p 

W = ^ 

Wt 

Expressing furthermore the speed of sound in the throat area by 
equation 98 finally yields 


h - hstat + 


Pt^ Vs, 


2 p 


tRtTt 

1 


(107) 

(108) 
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Applying this result to equation 98 yields the partial derivative 


shown below: 

^ ^ -/ ShstaA , PtSs. tRtTt iK?) I (109) 

According to equation 100 we obtain 

( ^^stat \ _ P 
\ainp;3 -p 


Furthermore 



Using the general gas equation, we obtain 



PtSs. tRT 1 

P^VsRt Tt e^ J 


( 110 ) 


Substituting the proper terms in equation 97 finally results in 


Pk+.=Pk 


1 + 


hch - (hstat + p 2 Tt S 




jZ -O T* . 9 ^ i 


)■ 


RT 


■ P2 RtTt\s ^ j 


( 111 ) 


MAXIMUM SPECIFIC IMPULSE 

Optimum rocket engine performance is achieved when the mixture 
ratio provides the maximum specific impulse. The theoretical perform- 
ance calculation in this report includes an option to determine optimum 
mixture ratio for a constant combustion chamber pressure. The 
following computation procedure is used. Initially, a mixture ratio 
is assumed, and the corresponding specific impulse is determined. The 
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mixture ratio is then increased twice by a chosen increment, and the 
pertinent specific impulse values are calculated. These t-hree points 
in an Isp - p. - diagram are the basis for the extrapolation routine 
explained below. 


The specific impulse as a function of mixture ratio for a constant 


chamber pressure can be represented by an equation 

I - F(f^) 
sp 
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Based on the known three points, a mean slope can be determined between 


two points 

• ^spi " ^sp2 
mi = — ^ ^ 

1^1 - [^2 


m. 


^spz ~ ^sp3 
^2 “ 


which also holds true for the mean mixture ratios 






^ 1^1 + 
2 


- fe + M >3 
2 


The terms m^ and m^ can be understood as the first derivative of the 
original function. A second derivative can be developed by taking into 
account the change of the slopes m 


mx 


mi " 
xi -X 2 


at a mean mixture ratio of 


y = 


xi + X 2 

2 


Since the specific impulse as a function of mixture ratio is represented 
by a curve that resembles an inverse parabola, the first derivative 
can be assumed to have a constant slope. Therefore, the equation of 
a straight line for the first derivative yields: 


m 


m 2 


mi - mz 

Xi - X2 


(x - X2) 
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Considering the second derivative, this equation changes to 


m - m2 = mi (x - X2) 

'The condition for the Tnaximutn specific impulse requires 


m = 0 


leaving 


- m2 = " ^2) 


The three originally calculated points are used to compute m2» > 

$nd Xj, A new approximation for the mixture ratio p. can be obtained 
by solving the latter equation for x. 


X 


miX2 - m 2 

mi 


( 112 ) 


At the point of maximum specific impulse, a secant changes into a tangent; 
the difference between two consecutive x-values approaches zero, or 

X ■ fJ. 

This approach exactly determines the maximum for a parabolic curve 
when three points are known. Since the Isp- P- curve is not an exact 
parabola, the predicted maximum is not identical with the real one. 
However, with every new calculation point a new set of three points 
is available to predict a new and better point. This approach physically 
moves an inverse parabola along the line Isp = F ( H-)p a constant 
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chamber pressure until its maximum converges into the real maximum. 

In the calculation procedure presented in Appendix G, it is assumed 
that a maximum specific impulse is obtained when 

< 0.01 

M-l - \^2 

To avoid any point of inversion, the mixture ratio determined for the 
maximum specific impulse is changed by Ap= 0.001 to the other side 
of the maximum to determine whether the specific impulse in this 
location is smaller than the previous one. 

THRUST LEVEL UPRATING 

When the thrust of an engine is uprated without changing its 
geometry, the new operating parameters, mixture ratio and chamber 
pressure, can be calculated in a separate option using the calculation 
program described in Appendix G. 

The thrust increase can be invoked by changing the chamber pres- 
sure without varying the mixture ratio. In this case a new chamber 
pressure can be determined under the assumption that the chamber 
pressure is directly proportional to the thrust 



^Des 


( 113 ) 
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The thrust value (F^) (calculated using the program) belonging to 

the new chamber pressure ) must be compared with the desired 

thrust (F ) . If both values do not fall within a certain 
Des 

extablished tolerance, a new chamber pressure ) will be determined 
by applying a straight line extrapolation. 


P - p + /p 

Ch2 Ch2 F2 - Fi ' 


Des 


-F2) 


(114) 


This procedure must be repeated until the desired thrust value 
is obtained (within a certain tolerance). 

Figure 11 indicates that lines for constant chamber pressure, 
representing the relationship between specific impulse and mixture 
ratio, show a shift in maximum specific impulses to higher mixture 
ratios with increasing chamber pressure. 
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FIGURE 11. SHIFT OF MAXIMUM SPECIFIC IMPULSE AS A FUNCTION 
OF CHAMBER PRESSURE AND MIXTURE RATIO 



For uprated engines a mixture ratio shift should also be taken into 
account to obtain high performance. 

One method is to determine the chamber pressure so that the desired 
thrust level is obtained, and the specific impulse is a maximum. This 
option is included in the computer program by combination of the maximum 
specific impulse and the constant mixture ratio uprating procedure. 

A second method considers that most rocket engines operate by some 
percent off the optimum mixture ratio, which yields maximum Isp. Since 
the mixture ratio of a rocket engine shifts during flight because of 
acceleration, propellant column height above pump inlet, density 
change of the propellants as a function of temperature, and many other 
reasons, it is preferred to stay on the side of the maximum impulse that 
shows the smaller slope. This allows a greater variation in mixture 
ratio for the same Isp tolerance, compared to the other side having a 
steeper slope. 

The calculation proceeds as follows: 

At first the combustion chamber at its original performance 
level is considered. Then for the same chamber pressure, the mixture 
ratio is determined that will yield maximum specific impulse. The 
difference between these mixture ratios is expressed as a percentage , 
and the sign indicates on which side of the optimum p. -value the 
combustion chamber is presently operating. Next the chamber pressure 
is varied, and the new point of maximum Isp, with the corresponding 
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optimum mixture ratio is calculated. Deviating from this optimum mix- 
ture ratio value by the initially calculated per cent value, an uprated 
mixture ratio is found, with which the desired thrust level and the cal- 
culated level must be compared. If the thrust level difference does not 
satisfy a chosen tolerance, the whole procedure must be repeated. 

COMBINATION OF SHIFTING AND FROZEN EQUILIBRIUM 
As stated in Reference 10, dissociation is becoming more and more 
pronounced when the temperature rises above 1500 °K. In the exit area 
of rocket engines, the static temperature is normally below this value, 

A real expansion process does not follow either the shifting or frozen 
equilibrium assumption, but occurs somewhere between these two cases. 

To calculate data that are close to real values, a combination of both 
concepts was considered in the calculation program and can be used if 
desired. The calculation starts with shifting equilibrium and switches 
to frozen equilibrium when the static temperature of the exhaust products 
falls below T = 1500°K. 

PARAMETERS 

Thermodynamic and rocket performance parameters are calculated 
according to the generally indicated procedures in text books (Ref. 4). 
The equations used in the computation process presented are indicated 
in Appendix G. 
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TOLERANCES 

When a combustion-expansion process is calculated, an inaccuracy 
of the results develops because of the influence of the two main factors. 
First,^the thermodynamic data, entropy, enthalpy, and specific heat at 
constant pressure, vary with different sources. Data for the heat of 
formation show discrepancies that are caused by the way they were 
obtained (calculated, measured, determined under the assumption of 
certain molecular structures) . In forming polynomials for these 
properties as a function of temperature, an additional inaccuracy deve- 
lops depending on the type of polynomial and the temperature range 
covered . 

Secondly, the calculation procedure does not produce exact answers. 
Numerical solutions and iteration procedures can only yield answers 
within a certain tolerance. This refers especially to the determination 
of the mole numbers, temperature, enthalpy, entropy and in particular 
to the throat and exit area conditions on which all the rocket perform- 
ance parameters are dependent. 

A theoretical error analysis for the complete calculation procedure 
is not favourable because of the complexity of the program. 

Usually the following tolerances resulted in a better accuracy 
than required (compare Reference 1). 

With the nomenclature of Appendix B we apply to: 

1. Equilibrium equations 

AF. < 5 X 10"^ 
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I-Ii 

a 

o 


1 - 




-7 

< 

5 X 

10 



-7 

< 

5 X 

10 


3. Pressure equation 

4. Enthalpy equation 

^ HO 

5. Entropy equation 

F 

1 - ^ 

so 


< 5 X 10 


< 5 X 10 


-7 


< 5 X 10 


-7 


etc. 


6, Tolerance for the throat area condition 

.-4 


1 

w 


< 2 X 10 


7. Tolerance for the exit condition 

.-3 


1 - 


' Des 


< 1 X 10 


These tolerances are too stringent in some cases and must be 
relaxed to obtain a solution. When a sequence of calculations is 
performed over a range of chamber pressures and mixture ratios, the 
tolerances must be relaxed so that the most critical case will pass. 

During the iterations solving for mole numbers and temperature, 
the calculated new approximations sometimes lead to a divergence and 
failure in the computations. This occurs occasionally when the first 
assumptions of the mole numbers were not close enough, and the 
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corrections caused divergence or an underflow in the electronic computer. 
To avoid such breakdowns, the following procedures were applied: 

1. If the mole numbers are corrected and become so small 

that the computer assigns them equal to zero, these 

-35 

terms shall be reset to a small value (1 x 10 ). 

2. If .the corrections become greater than 10, they are to 
be divided into 10. All corrections shall then be 
multiplied by the maximum absolute value obtained. 
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APPENDIX A 


THERMODYNAMIC RELATIONS 

Each property of state can be expressed by two other properties 
of state. From the relationship F (P, v, T) = 0 an important equation 
is derived. 


V = 

T = 
dv = 


dT 


F (P, T) 
F (P, v) 


(-1^) - ^ (i?-) 


dP 


dv 


( 116 ) 

(117) 


V P 

Replacing dv in one of the equations and rearranging the terms yields 

-1 = 1-^1 1-^1 (118) 


■ (4t), m. Hi) 


A similar expression will be developed from the relationships 

T - F (u. V) 
u = F (T, v) 


dT 


du = 


= (-4^) ^ 

- (■#) - * (t7) 


dv 


dv 


(119) 

( 120 ) 

( 121 ) 


( 122 ) 


Substituting du in one of the equations and rearranging gives 

1 = (-^1 (123) 


m. m. m 
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The first law of thermodynamics is represented by the equation 


dq 

= du + Pdv 

(124) 

The second law of 

thermodynamics reads 


dq 

= Tds 

(125) 

Introducing the definition of enthalpy 


h 

= u + Pv 

(126) 

Helmholtz function 


a 

= u - Ts 

(127) 

Gibbs function 



g 

= h - Ts 

(128) 

their differentials 

are 


dh 

= du + Pdv -1- vdP 

(129) 

da 

= du - Tds - sdT 

(130) 

dg 

= dh - Tds - sdT 

(131) 

Application of the corresponding equations 
second law of thermodynamics yields 

representing the first and 

dh 

= Tds -h Pdv 

(132) 

da 

= -Pdv - sdT 

(133) 

dg 

= vdP - sdT 

(134) 


Futhermore, the following relations hold true 
a = F (v, T) 
g = F (P, T) 
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u = F (s, v) 

h = F (s, P) 

s = F(T, P) 
s = F (T, v) 

T = F (v, P) (135) 

Their derivatives are 


da = 

Or) 

dv + ( 
T 

dT 

, 3T/ 

V 

(136) 

dg = 

m) 

1 dP + 
T 

(ifl 

(137) 

du = 

Us; 

1 + ( 
V 

■if) - 

S 

(138) 

dh = 

(4^) 

ds + 1 
P 

dP 

^ a p/ 

s 

(139) 

ds = 

/as' 

V a T . 

) dT + 
P 

(iw), - 

(140) 

ds = 

1 d s \ 

1 a T > 

1 dT + 

V 


(141) 

dT = 

/_a_T_' 

\ a V y 

) dv+ I 
p 

(44) - 

V 

(142) 
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Comparing the equivalent terms in the corresponding equations 


(a:) 

1 = T 

V 

(113) 

(4t) 

1 = 
s 

(144) 

(-^1 

1 = T 

(145) 

\ d s 1 

f 

P 


(-U) 

1 = V 

T 

(146) 


1 = "S 

(147) 

\ 

P 

/ah' 
\ a p > 

1 = V 

s 

(148) 

/ a a 
\ a T ^ 

1 = 

(149) 


All of the above equations involve properties and are exact differentials. 
Since properties are independent of a path, the following mathematical 
concept applies. 

dz = Adx + Bdy 



This results in the subsequent relationships 



s V 


( 150 ) 
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( 151 ) 



The specific heats are, by definition 


Cv 

/ 9u 

1 ax 

1 

Cp 

/ 9 h 
\ a T 

•)p 

Expanding these two equations indicates 

c 

V 

/ au 

\ a s 

■) Ht) 

V V 

Cp 

/ 9h 

\ a s 

■) (-#) 

p p 

Applying equations 

(143) and 145) and expanding the terms 

Cv 

= T (- 

as \ T ( W ap\ 

9T' \ ap / ' a T / 

V V V 

Cp 

“ ^ (- 

as\ Y /3s\ ( \ 

9T ^P ^ 9v /p\ ax /p 

Substitution of terms from 

equation (150) and (151) yields 

/ ap 

- ■ 

/ ax\ 

V 9T 

/ 

V 

X \dvf^ 

/ av 

- 

^p ( ax \ 

V aT 


X 1 ap /g 


(152) 

(153) 

(154) 

(155) 


(156) 

(157) 


75 



APPENDIX B 


APPLICATION OF NEWTON- RAPHSON METHOD TO 
SOLVE FOR MOLE NUMBERS AND TEMPERATURE 

This appendix applies the Newton-Raphson method to equations 
39» 40, 41, 42, 44, 45, and 46 to linearize the original non-linear 
equations that are used to solve for the mole numbers, N^, and the 
temperature, T. 



1 


Equilibrium equation 


F. = 0 = 



+ 


In Pi 


- ailnP^ 


- bilnPy - CilnPjj 


(158) 


AFi = 0 - Fi for the right solution Fi = 0, otherwise Fi has a 

value 


8Fi 8Fi 9Fi 

AFi = -tt- AA + APi + AT 
^ 8 A 8 Pi ^ 8 T 



8 Fi 

1 

8 Fi _ 

ai 

1 

•H 

CO 

bi 8 Fi 

Ci 

8 Pi' 


8 P 2 

■ Pz’ 

9 Py 

■Py’ aPx 

" ■ Px 

8 Fi 

8 

/agt°\ 

8 

/aht° - 

ATS°\ / 

AHt°\ 

8 T 

’ 8 T 

Vrt /i 

" 8 T 

\ RT 

H' 

RT^ ) 


1 

a-i 


bi 

•=i . 


AFi: 


^i — 

APz 

AP, 

ir'y 

7 - ^ 

**•••• 



AFi = A In Pi - aiAlnP 2 - biAlnPy - CiAlnP^ - 

/aht°\ 

- I — I AlnT 

\rt 


(159) 


When a condensed product is formed, the partial pressure 
must be neglected due to the section entitled "condensed phases. " 
The equilibrium equation for condensed phases yields: 


AFi = “ aiAlnPz - biAlnPy - CiAlnPx 



(160) 


2. Mass balance equation for element Z 


1 

^a = a-tot - Xrf ^i^i 
^ i=l 


(161) 
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^-tot - 


for the correct solution 


1 

~ ~T ^2 
^ i= 1 


9f'a 


AA + — ANi + — AT 


9Fa 

TT 


1 

I' 


115 

L aiNi 
i= 1 


9Fa _ 9Fa 
9Pi “ 8Ni 



9 Fa 
9T 


= 0 


AF = 
a 



1 1 

^Fa = - ;^ E aiNiAlnA + — ^ aiNiAlnNi 
i = 1 i = 1 


( 162 ) 


3. Pressure equation 


= P = 


110 


F Pi 

i= 1 


P = Pq for the correct solution 


no 

AFp = Po - E Pi 


(163) 
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9Fp> 9 F p) 9Fp 

^^P= — 


9Fp _ 

^ i= 1 ' 


no no 

AFp = Y. 1 = L Pi^lnPi (164) 

i= 1 i= 1 

The summation of partial pressures can only be performed 
for the gaseous products because the partial pressures for condensed 
phases have to be neglected. 

4. Enthalpy equation 

115 

Fh=H=xE (HtVi (165) 


H = Hq for the correct solution 


X ID 

AFtt= Hq - X E (Pt ^i^i 


9Fh 9PH 9Fh 

+ — AT 


X 4. 

H 1 V - 0. 
r = - T7 E (Ht )iNi 


9A A 
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9 Fh 9 Fh 1 .Q 0. 

2^ (Ht )i 


8P 8N A . 


i= 1 


115 


115/ \ 115 

i ‘f (HT°)iNi = A E^(/cp”dTjjNi = 


8T 8T / , 
1=1 


( MNi 


115 


115 


AFpj = - ^ Z] A 


i= 1 


i= 1 


115 


_ 0 


+ E (Cp)iNjAT 
i= 1 


, 1*5 , 115 

iI'Fh^-xE (HT')iNi AlnA + i X (Hx“)i Nj A In N; 


i = 1 


i= 1 


115 


+ T L (Cp®)iNiAlnT 
i= 1 


5. Entropy equation for gaseous products 


no 


Fc = S =x E [(ST°)i - R InPilNi 


i= 1 


S = SO for the correct solution 


no 


AFc5 = so - E ^ [(St°)i - R In Pi] Ni 


i= 1 


8Fo 8Fs 9Fs 

^^S=-^^A + — APi + — AT 


( 166 ) 


(167) 


80 


E [(ST°)i - RlnPilNi 


3Fs 8Fs 1 "0 ,, 0 , 1 /"» 


^ = ^ = A E 'V)i ■ il E SlnPi+ N, E 
^ ^ 1=1 \l=l 1=1 


8Fs 8 
8T ■ ‘8T 




no no 

= - TI E [(St°)i - RlnPjNi AA + - ^ (ST°)i ANi 


A" ,r 


-i: RlnPiANi - E Nil: ^Pi + X E Ni(-|-). AT 

i = l i=l ^ i=l\/^ 

no fiio 

AFs= -i E [(ST°)i - RlnPi]NiAlnA + - E (ST°)i Ni A In Ni 
i = 1 i = 1 

no no _ 

- E RNilnPi AlnNi - E RNiAbiPi 

i=l i=l 

1 

+ T E Ni(Cp°)i AlnT (168) 


For condensed phases, the partial pressures must be 
neglected. The entropy equation is then 


1 115 

S = A E (ST)iNi 
i= 111 


(169) 
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Applying the Newton-Raphson method, we finally obtain; 


115 115 

AFs=x E N.(Cp°)iAlnT + - E Ni(ST°)j 
i= 111 i= 111 


1 

-T E 


(ST°)iNi AlnA 


4=111 


Ain 


( 170 ) 



APPENDIX C 


APPLICATION OF NE WTON-RAPHSON METHOD TO SOLVE 
PARTIAL DERIVATIVES FOR 
CONSTANT PRESSURE AND TEMPERATURE 

This appendix appli-es the Newton-Raphson method to equations 39» 

40, 41, 42, 44, 45 and 46 to linearize them so that the partial derivatives 
for constant pressure and temperature can be solved easily. 
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Derivatives at Constant Pressure 


1. Equilibrium equations considering gaseous compounds expressed 
by partial derivatives at constant pressure: 



0 


AHt \ , d In Pi _ dlnPz ^ dlnPy ^ d lii Px 

"f^l. dlnT ■ dlnT " dlnT "^^dlnT 

(172) 


2. For condensed phases the equation is identical except that the 
term 


din Pi 
dlnT 


= 0 


0 = 



dlnPz , dlnPy dln Px 

dlnT " dlnT " dlnT 


(173) 


3. Mass balance equation expressed by partial derivatives at 
constant pressure ; 



i = l 


(174) 
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115 

0 = -aQA+ ^ ^ 

i = 1 


0 


-ao dA 



aj dN^ + 



dao and da^ cannot change since mass is constant. 


115 


0 = -aQ A d In A ... . + a j d In Nj 

i=l 


din A . .. dlnNj 

dlnT 


- . dlnA . 

0 = -a.Q -p • • • ♦ + ^ ^i ^i 


i = l 


or 


115 115 

0 = - ajNi dlnA + ^2 ^i^i dlnNi 
i=l iVl 


115 ,, . 115 

“ = -.E E 

1=1 i = 1 


dlnNi 

dlnT 


(175) 


4. Total pressure equation expressed by partial derivatives at 
constant pressure: 


P = 


110 


E Pi 

i - 1 


(176) 


110 

8 E Pi 

9P 8Pi 


Applying the mathematical 


dPi 

concept 


dlnx _ J. 
dx X 
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dx = X d In X 


we obtain 

no 

Pi d In Pi = P d In P 

iVl 

at constant total pressure the d In P term is equal to zero. 
P d In P = 0 


Dividing the equation by d In T it follows 


0 = 


no 


i: Pi 

i - 1 


din Pi 
dlnT 


Derivatives at Constant Temperature 


5.. Equilibrium equations expressed by partial derivatives at 
constant temperature^ 

/agT\ 

0 = 1 —: 2 :: I + In Pi - ai In P^ - bi In Py - Ci In P^^ - . • . . 

' ' 1 ( 


0 = d 


AGt 


+ d In Pi - aidlnP 2 - bidlnPy - CidlnP^ 


AGt 


/i d InPj d InPz . . d In Py d InP, 

dlnA ^ dlnA dlnA ^ d InA d InA 


dlnA 


is only a function of temperature 



d InPz 
d In A 


d In Px 
d InA 


(179) 


0 = 


d InPj 
d InA 


d InP 


y 


^ d InA 


6' Mass balance equations expressed by partial derivatives at 
constant temperature; 


115 

^o ~ y ^i ^i 

i= 1 

115 

0 = -Aao + 2 
i= 1 


(180) 


Since = Pj^ we obtain 


A 


8aQ 

9 aQ 


dao+ ao 


9A 

9A 


dA 


115 




9Nj 

9Ni 


115 

dPi +^Ni 
i- 1 


9aj 

8ai 


da ^ • 


Because of the constancy of elements 


dao = 0 
da^ = 0 • 

Dividing the equation by din A and applying the following mathematical 
concept 


d Inx ^ 2. 
dx X 


or 


dx = X d In X 


the result is 


0 = -a^A 

0 = -UqA 


115 


d InA 
d InA 

i = 1 
115 

+ VaiNi 


N, 


d In Pi 
1 d InA 


d InP 


i = 1 


I 


d InA 


(181) 
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APPENDIX D 


DERIVATION OF THE CHANGE IN MOLECULAR 
WEIGHT AS A FUNCTION OF PRESSURE AND 
TEMPERATURE 

Performing a mass balance due to the reaction equation. 

AlNj.(Z^p, Ybp, XcF '’be - ^c6 >1 

= SNi( Z^i, Ybi. 
we find for each element 

A(Npap + N0 ag) L^ = ^N^a^L^ 

A(Npbp + Ng bg) Ly = EN^b^Ly 
A(NpCp + Ng Cg) L^ = SNm.L^ 

A[Np(apLz+ bpLy+ CpL^^ + •.••)+ Ng (ag L^ + bg Ly + CgL^ 

+ ....)] == 2N^(a^Lz f" bj^Ly + c^L^ + 

In general the molecular weight is defined as 

M. = a.L„ + b.L„ + c.L„ + (184) 

1 1 Z 1 Y 1 X 

In a reaction process, the mole numbers and the molecular 
weights of the reactants are known. This means the relative mass is 
a constant. The mass balance equation reads now 


1 15 



i= 1 


Very often the mole numbers of the reactants are determined so 
that the total mass of the reactants is equal to one. In this case, 

mj. = 1 . 


(182) 


(183) 
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An average molecular weight can be defined for a mixture as 

115 

NM = Yi NiMi 
i= 1 

115 

Y NiMi 

. , i= 1 


no no 

With N=^ Pj = P (neglecting the volume of the 

i = 1 i = 1 

condensed products) and the total mass of the products equal to the 
total mass of the reactants, we obtain 


M = 


A mj. 
P 


( 185 ) 


Rewriting this equation in a logarithmic form and differentiating 
with respect to In T results in 

dlnM _ din A dlnmp dlnP 
d In T ” dlnT d lnT “ dlnT * 


Since the mass, mj!, will not change. 


d In mj. 
dlnT 


= 0 . 


For a constant pressure process, P does not change and 
dlnP ^ 


Finally we obtain 

( dlnNt\ _/dlnA\ 
dlnT Jp “ I dlnxlp ‘ 
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starting again from equation (185) in logarithmic form but now 
differentiating with respect to InP yields 

d In M _ d In A d In m r d In P 
dlnP ~ dlnP d In P ‘ dlnP * 

Due to the constancy of mass 


d In m r _ 
dlnP “ 

The final result is 

d In M _ d In A 
dlnP “ dlnP ‘ 

and can also be applied to a constant temperature process 
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APPENDIX E 


CONSTANTS AND 
1.98726 

2.7182818285 

32.1740 

3087.16928168 

2.2046 

14.696 


CONVERSION FACTORS 


Universal gas constant 


Euler number 


Gravitational constant at sea level 


ft/lb 


lb 


lb 
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APPENDIX F 


A 

ACC 

AISP 

AISPX ) 
AISPY > 
AISPZ ) 

AISPM 

AISPRA 

AISPR-V 

ALPHA 

AMNF 

AMNFS 

AMN0 

AMN6S 

AMNRi 

AMRW 

AMRWD 


NOMENCLATURE FOR CALCULATION PROGRAM 


Multiplier in reaction equation to make tht 
partial pressures and mole numbers numerically 
equal 

Injector plate area 


Local specific impulse 

Previous values of AISP during iteration for 
maximum sped f ic impulse at constant chamber 
pressure 


sec 

sec 

sec 


mole 

mole 

[mole 


Maximum specific impulse 

Real specific impulse at predetermined altitude 
Real vacuum specific impulse 

Mole per cent value of oxidizer A in a mixture 

Mole number of the fuel, determined as one (1) 

Mole number of the fuel to yield a molecular 
weight of one (1) for the propellants 

Mole number of the oxidizer to satisfy the required 
mixture ratio of the propellants by weight 


mole 


1 


Mole number of the oxidizer to yield a molecular 
weight of one (1) for the propellants, considering 
also the desired mixture ratio by weight 

Mole ratios of the combustion products 


Mixture ratio by weight 

Chosen mixture ratio increment for iteration on 
maximum specific impulse at constant chamber 
pressure 
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AMRWL 


AMRWM 


AMRWS 

AMRWX ) 
AMRWY > 
AMRWZ ) 

AMW 


AMWP 


AMWR 


ANA 

ANB 

AND 

ANG 

ANi 

AS 

ATWC '' 

ATWF 

ATWH 

ATWHE 

ATWN ^ 

ATWe 

ATWU 

ATWX 

ATWY 

ATWZ ^ 

BETA 

Ci 


C* 


[-] 

H 

H 

H 


Original mixture ratio by weight to determine the 
PERCT value for thrust level uprating 

Optimum mixture ratio by weight yielding maximum 
specific impulse for considered chamber pressure 

Stoichiometric mixture ratio 

Previous values of AMRW during iteration for 
maximum specific impulse at constant chamber 
pressure 



Average molecular weight per mole of gaseous 
products 

Relative mass of propellants for the required 
mixture ratio 



Relative mass of the reaction products defined as 
one (1) 



Number of Nitrogen gram atoms in oxidizer A and B, 
and Fuel G and D 


Number of gram atoms in the combustion products 
Throat area 



Atomic weight of Carbon, Fluorine, Hydrogen, 
Helium, Nitrogen, Oxygen, with U, X, Y, Z as 
available spaces for other elements 



cent value of oxrdrzer B rn a mrxture 

Number of carbon gram atoms in the combustion 
products 

Characteristic velocity 


93 



CA 

CB 

CD 

CG 

CF 

CFA 

CFRA 

CFRV 

CFV 

CKl 

CK2 

CK3 

CK4 

CN^ 

CNS 

CNSMCP 

CONVER 



mole 

[mole 



[-] 


CPP 

C*R 

CRS 

CRSS 

cvv 

DELTA 



cal I 
mole°K 

W] ■ 



Number of carbon gram atoms in oxidizer A and B, 
and fuel G and D 

Local thrust coefficient 

Thrust coefficient at predetermined altii 

Real thrust coefficient at predetermined altitude 

Real vacuum thrust coefficient 

Vacuum thrust coefficient 

Conversion factor 

Conversion factor 

Conversion factor 

Conversion factor 

Mole numbers of the combustion products 

Summation of the mole numbers of the combustion 
products 

Summation of the mole numbers of the gaseous 
combustion products only 

Error left in the iteration process due to the 
approximated values 

Specific heat at constant pressure of the individual 
combustion products 

Specific heat at constant pressure considering all 
combustion products 

Real characteristic velocity 

Local isentropic exponent 

Isentropic exponent in the throat area 

Specific heat at constant volume considering all 
combustion products 

Mole per cent value of fuel D in a mixture 
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DlnA 

DlnP^ 

DPA 

DPM 

DPP^ 

DQi 

DIM 

DTPi 

E 

EQUIVR 

ETAC* 

ETACFV 

ETAW 

EXPR 

EXPRL 

Fi 


H 

H 

H 

H 

H 


[-] 

H 

[-] 

[-] 

[-] 



Enthalpy difference between the combustion products 
and the elements from which they are formed. All 
values divided by R*T 

Finite difference of the logarithm of the reaction 
equation multiplier. 


Finite differences of the natural logarithm of the 
partial pressures (mole mumbers) 

Change of molecular weight as a function of 
temperature* at constant pressure 

Change of molecular weight as a function of 
temperature at constant pressure 


Partial derivatives 



at constant pressure 


Constant term in the equilibrium equation 


Change of molecular weight as a function of pres- 
sure at constant temperature 


Partial derviatives 
temperature 



at constant 


Euler number (2.7182818285) 


Equivalence ratio; mixture ratio divided by the 
stoichiometric mixture ratio 

C'orrection coefficient for the characteristic 
velocity 

Correction coefficient for the vacuum thrust 
coefficient 


Correction coefficient for the mass flow rate. 
Expansion ratio 
Required expansion ratio 

Number of fluorine gram atoms in the combustion 
products 


95 



FA 



FB 


g 

FD 

r 

nole 

FG 


FL 

ibj 

FROZ 

— 

FT 

lb 

FTA 

lb 


k « 

FTRA 

lb 

B X 

FTRV 

lb 

FTV 

lb 

FX ) 


[ihl 

FY f 


rJ 

GAMMA 

— 

H 

cal 

Hi 

-1 

mold 

HA 


HB 



HD 

HG 


mo 1 ej 

HEi 

‘ g 

mole 

HEA 



HEB 


‘ g 

RED 


mole 

REG 


RO 

f] 

RTi 

mole 

RTFD 

1 

cal ' 

RTFG 

f 

mole^ 


Number of fluorine gram atoms in oxidizer A and 
B and fuel D and G 

Desired thrust level 

Switch constant for frozen equilibrium 
Local theoretical thrust 

Theoretical thrust at predetermined altitude 
Real thrust at predetermined altitude 
Real vacuum thrust 
Theoretical vacuum thrust 

Previous values for FTRA during iteration for 
desired thrust level 

Mole per cent value of fuel G in a mixture 

Enthalpy of all combustion products 

Number of hydrogen gram atoms in the combustion 
products 

Number of hydrogen gram atoms in oxidizer A and B, 
and fuel D and G 

Number of heliuin gram atoms in the combustion 
products 

Number of helium gram atoms in oxidizer A and B, 
and fuel D and G 

Enthalpy of the propellants 

Molar enthalpy of the combustion products 

Enthalpy of fuel D and G at normal boiling point 
or at T = 298.15^ K 
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HT0A ) 
HT0B f 

Ik 

IB 

IC 

ID 

IE 

IF 

IG 

IH 

II 

IJ 

IK 

IL 

IM 

K 

KE 

KF 

K0 


"cal ^ 
mole 

Enthalpy of oxidizer A and B 
point or at T = 298.15° K 

at normal boiling 

— 

Switch 

constant 

for throat area calculation 


— 

Switch 

constant 

within the throat area calculation 

— 

Switch 

constant 

considering 

either 

calculation at 


a requi 
downs tr 

red expansion ratio or many 
earn of the chamber 

local conditions 

— 

Switch 

constant 

considering 

static 

pressures 

in a 

L. _i 

chosen 

range 





— 

Switch 

constant 

considering 

static 

pressures 

in a 

- 

chosen 

range 





— 

Switch 

constant 

considering 

static 

pressures 

in a 

- 

chosen 

range 







Switch 

cons tant 

considering 

the calculation 

of a 


desired thrust 

level 




1 

Switch 

constant 

considering 

the calculation 

of a 


ft 


sec_ 

mole 


mole 


maximum specific impulse 

Switch constant within the calculation routine 
for maximum specific impulse 

Switch constant within the calculation routine 
for maximum specific impulse 

Switch constant considering the difference between 
an optimum and required mixture ratio 

Switch constant considering a combination of a 
desired thrust level and a maximum specific impulse 

Switch constant within the desired thrust level 
calculation routine 

Test constants 

Gravitational constant at sea level, 32.174 
Molecular weight of fuel 
Molecular weight of oxidizer 
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KKli ^ 

KK2i 

KK3i 

KK4t ► 

KK5i 

KK6i 

KK7i 

6i 



mo 


le_ 


Polynomial coefficients for calculation of 
enthalpy, entropy, and specific heat at constant 
pressure 


Number of oxygen gram atoms in the combustion 
products 


0A 

6B 

6D 

0G 

P 


PC 

PAM 

PCM 


g 

mole_ 

atm 
L = 

atm 

atm 

atm 

atm 


PERCT 


PO 

PS 


atm 


atm 


Number of oxygen gram atoms in oxidizer A and B , 
and fuel D and G 


Local static pressure of all combustion products 

Local partial pressures of the individual 
combustion products 

combustion chamber end stagnation pressure 

Ambierut pressure at considered altitude 

Optimum combustion chamber end stagnation pressure 
which belongs to a required thrust level yielding 
maximum specific impulse 

Percent value considering the deviation in mixture 
ratio off the optimum value yielding maximum 
specific impulse 

Predetermined local static pressure of the 
combustion products 

Static pressure of the combustion products in the 
throat area 


PSUM 

PX ) 
PY / 

R 

RS 


atm 

atm 

cal 

mole 

cal 

_g-°K 


Local static pressure of all combustion products 

Previous values for PC during iteration on desired 
thrust level 

Universal gas cons-tant 1.98726 

Specific gas constant for all combustion products 
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RSS 


cal 

iT^K 


Specific gas constant for all combustion 
products in the throat area- 


c 

cal 


D 

g°K 


SHIFT 



r 

SO 

cal 


ST^ 


■ cal 


mole 

T 




_ 


TL 



TC 



TS 

°K 


Tang a 

_ _ 


Ui 


g 

-| 


mol< 


UA^ 


UB 


r ^ 

“ 

UD 


mol* 

e 

UG . 


V 

ft ■ 
sec 


VC ^ 


VF 


VH 


VHE 


VN 




ve 

1 



vu 


vx 


VY 


vz ^ 


V3 


1 


sec 


Entropy of all combustion products 

Switch constant considering shifting equilibrium 
calculation 

Calculated entropy value which is considered in 
the isentropic expansion process 

Molar entropy of the combustion products 

Temperature 

Temperature limit counsidering the transition from 
shifting to frozen equilibrium calculation procedure 

Theoretical combustion chamber temperature 

Theoretical static temperature in the throat area 

Slope of the secant for maximum specific impulse 
calculation. 

Number of gram atoms for an element (to be 
determined) in the combustion products 


Number of gram atoms for an element (to be 
determined) in oxidizer A and B, and fuel D and G 


Local velocity 

Most common valence of carbon 

Fluorine 

Hydrogen 

Helium 

Nitrogen 

Oxygen 

Element to be determined 
Element to be determined 
Element to be determined 
Element to be determined 

Local speed of sound 
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vss 

ft I 


sec 

w 

' lb ' 


sec 

WR 

' lb " 


sec 

Xi 

r “I 

g 


mole. 

XA ^ 


XB 


g 1 

XD 


mole 

XG , 


Yi 

g 1 


mole 

YA ' 



YB 

► 

§ 1 

YD 


mole_ 

YG . 



Zi 

ZA 
ZB 
ZD 
ZG , 


g 

mole 

' g 
mole 


OC 

OF 

OH 

OHE 

ON 

oe 

ou 

ox 

OY 

oz 


1C 

IF 

IH 

IHE 

IN 

le 

lu 

IX 

lY 

IZ 




> 


J 




speed of sound in the throat area 
Theoretical weight flow rate 
Real weight flow rate 

Number of gram atoms for an element (to be 
considered) in the combustion products 

Number of gram atoms for an element (to be 
considered) in oxidizer A and B and fuel D and G 

Number of gram atoms for an element (to be 
considered) in the combustion products 

Number of gram atoms for an element (to be 
considered) in oxidizer A and B and fuel D and G 


Number of gram atoms for an element (to be 
considered) in the combustion products 

Number of gram atoms of an element (to be 
considered) in oxidizer A and B and fuel D and G 


Number of gram atoms yielding a relative mass of 
one (1) for the propellants considering carbon, 
fluorine, hydrogen, helium, nitrogen, oxygen, and 
elements to be determined (U, X, Y, Z) 


Number of gram atoms yielding a relative mass of 
one (1) for the combustion products, considering 
carbon, fluorine, hydrogen, helium, nitrogen, 
oxygen and elements to be determined (U, X, Y, Z) 
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APPENDIX G 


CALCULATION PROGRAM 

In this section the complete calculation procedure is documented. 
Input data, determination of chemical elements and compounds, equa- 
tions, tests, option to solve for various desired conditions and output 
data are listed. The simplified flow chart in Appendix H shows the 
arrangement of all equations and tests. 

The computer listing and pertinent flow chai :s are available 
and can be obtained on request from the author. 
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GENERAL COMBUSTION PROGRAM 


1. Molecular Weight of Oxidizer [g/mole] 

2. Molecular Weight of Fuel [g/mole] 

3. Mole Number of Oxidizer to Satisfy Desired Mixture Ratio [mole] 

4. Relative Weight of Propellants for the Required Mixture Ratio [ g] 

5. Number of Gram Atoms Yielding a Relative Weight of One for the 
Propellants [ ~ ] 
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K0 = ALPHA(CA • ATWC + FA • ATWF + HA • ATWH + HEA • ATWHE 

+ ANA • ATWN + 0A • ATW 0 + UA • ATWU + XA • ATWX + YA • ATWY 
+ ZA • ATWZ) + BETA(CB - ATWC + FB - ATWF + HB - ATWH 
+ HEB • ATWTHE + ANB • ATWN + 0B • ATW0 + UB • ATWU 
+ XB • ATWX + YB • ATWY + ZB • ATWZ) 

KF = GAMMA(CG - ATWC + FG - ATWF + HG - ATWH + HEG - ATWHE 

-I- ANG • ATWN + 0G • ATW 0 + UG • ATWU -f XG • ATWX + YG ■ ATWY 
+ ZG • ATWZ) + DELTA(CD - ATWC + FD • ATWF + HD • ATWH 
+ HED • ATWHE + AND • ATWN + 0D • ATW0 + UD • ATWU + XD • ATWX 
+ YD - ATWY + ZD - ATWZ) 

AMRW-AMNF-KF 

AMN8 

AMWP = AMNF • KF + AMN0 • K0 

or. - AMN0(ALPHA- CA -f beta • CB) + AMNF(GA MMA • CG + DELTA • CD) 

AMWP 

OF _ AMN0(ALPHA • FA + BETA • FB) + AMNF(GAMMA • FG + DELTA • FD) 

AMWP ; 

OH - AMN0(ALPHA - ha -f BETA • HB) -f AMNF (GAMMA • HG + DELTA • HD) 

AMWP 

Ohe = AMN0(ALPHA- HEA -I- BETA - HEB) + AMNF(GAMMA- HEG + DELTA - HED) 

AMWP ^ 

OH _ AMN0(ALPHA • ANA + BETA - ANB) + AMNF(GAMMA - ANG + DELTA - AND) 

AMWP 

O0 - AMN0(ALPHA • 0A + BETA - 0B) -t- AMNF (GAMMA • 0G + DELTA - 0D) 

AMWP ^ 

OU - ^^^6(ALPHA - UA + BETA - UB) + AMNF(GAMMA • UG + DELTA - UD) 

AMWP ~~ 

_ AMN0(ALPHA - XA + BETA - XB) -t- AMNF(GAMMA • XG + DELTA • XD) 

AMWP 
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6. Mole Number of Oxidizer to Yield Molecular Weight of O- ^ jr 
Propellants [mole/g] 

7. Mole Number of Fuel to Yield Molecular Weight of One for Propellants 
[mole/g] 

8. Stochiometric Mixture Ratio [ - ] 


9. Equivalence Ratio [ - ] 

10. Enthalpy of Propellants [cal/g] 

11. Test on T for Proper Selection of Temperature Ranges for 
Entropy, Enthalpy and Specific Heat Polynomials. 

12. Local Enthalpy of the Combustion Compounds [cal /mole] 
i = 1 - 115 


13. Local Specific Heat at Constant Pressure of the Combustion Compounds 

[ cal/mole ®K] i = 1 - 115 

14. Local Entropy of the Combustion Compounds [cal/mole“K] 
i = 1 - 115 

15. Test Whether Shifting or Frozen Equilibrium is Invoked 
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AMN6(ALPHA • YA + BETA • YB) + AMNF(GAMMA • YG + DELTA ■ YD) 
" AMWP 

_ AMN6(ALPHA • ZA + BETA • ZB) + AMNF(GAMMA • ZG + DELTA • ZD) 
" AMWP 


AMN0S = 


AMN0 

AMWP 


AMNFS = 


AMNF 

AMWP 


AMRWS = (GAMMA [CG - VC + FG ■ VF + HG • VH + HEG -VHE + ANG • VN 
+ 0G -V0 + UG -VU + XG VX + YG -VY + ZG -VZ] 

+ DELTA [ CD • VC + FD • VF + HD • VH + HED • VHE + AND • VN 
+ 0D -V0 + UD -VU + XD -VX + YD -VY + ZD -VZ]) • K0/ 
-(ALPHA [ CA • VC + FA • VF + HA • VH + HEA • VHE + ANA • VN 
+ 0A-V0 + UA- VU + XA-VX + YA -VY + ZA *VZ] 

+ BETA [CB - VC + FB - VF + HB - VH + HEB -VHE + ANB • VN 
+ 0B • V0 + UB ■ VU + XB • VX + YB • VY + ZB • VZ] ) • KF 


EQUIVR = 


AMRW 

AMRWS 


HO = AMN0S(HT0A • ALPHA + HT0B • BETA) + AMNFS (HTFG • GAMMA 
+ HTFD> DELTA) 


HTi = 


RT 


•^KKli + KK2i + KK3i ^ + KK4i ^ + KK5i ^ + KK6i 



CPi = R-(KKli + KK2iT + KK3i T^ + KK4i T^ + KK5i T^) 


rnC. rpj rp4 

STi = R(KKli In T + KK2i T + KK3i j- + KK4i — + KK5i — + KK?i) 

If FROZ = 1 Skip to Step 20 

If SHIFT = 1 Continue 
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16. Simultaneous Set of Equations Solving for A, T, Pi 

a. Equilibrium Equations [ - ] 

Gaseous Products: i = 1 to 100 

Condensed Products: i = 111 to 115, in this case In Pi = 0 
The term of the left-hand side of this equation is equal to -DQi 
and approaches zero for a solution. 

The positive parenthetical term on the right hand side of the 
equilibrium equation, not multiplied by D In T, is abbreviated 
by DHi 
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-fs), -m..- ■ “-(i?-),.. • “• (S'), 

- AN. 


04 


^‘(^)uo : (f )i " " =i[(f )... - " ^{(f ).« ■ 


^ "{(f) 


- lnP103 


103 


+ e 


In PI 06 


+ Y 




{(f ).« - 


] ^ 

] ^ "{(f 




In PI 04 


- lnP107 


j + AN; 

Mm 


■im 


- lnP105 


105 


] 


108 


- InPlOS 


]"^{(f)n.-'"""°]} = 


Din Pi - CiDlnPlOl 


- FiDlnPlOZ - Hi'DlnPlOS - HEiDlnP104 - ANi'DlnP105 - 0iDlnPlO6 

- UiDlnPlO? - Xi'DlnPlOS - Yi'DlnP109 - Zi'DlnPllO 


^104 


“■m 

.(s),.j 


105 


DlnT 
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Mass Balance Equations [ 


] 


Total Pressure Equation [atm or mole] 

The partial pressures and mole numbers are equivalent with the 
application of A. 
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d. Enthalpy Equation [ cal] 

Replace enthalpy by entropy equation (l6e) when step 44 has been 
passed the first time. 


e. Entropy Equation [cal/“K] 

Use step l6d during the first calculation until step 44 is passed. 


17. Test on Tolerances 


During the first calculation loop 1 


S 

SO 


must be by-passed. 


During all consecutive calculation loops 1 


H 

HO 


must be by-passed. 


P, H, S, IC, IF, IH, IHE, IN, IG , lU, IX, lY, IZ, are the constants 
calculated during the iterations in the simultaneous set of equations 
solving for A, T, and P| (see step l6b, l6c, l6d, l6e), which have 
approximately the same value as PO, HO, SO, OC, OF, OH, OHE, 
ON, 06, OU, OX, OY, OZ. 
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K =- R . T ^ . Pi 


D In A 




R • T • D In Pi + R • T 


1 ("). '■] 


DlnT 


'i i [® i • Pi - ^i] - i • « : T |(f ) i • Pi] } 

+ Rp£i(^)-Pi DlnT - Pi(l + InPi) Pij-DlnPi 

■ DlnPi-R-J^ (f),.Pi-Pi.lnPi].DlnA 


1 - ^ < 5 X 10-’ 




1 - ^ |<5xl0-’ 
‘ ^ < 5 X 10-’ 


1 .^|<5xl0-’ 

1 - ^ < 5 * 10-’ 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 


Continue otherwise go to step 18 
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DQ^ is the constant in the equilibrium equation which approaches zero. 


18. New Estimates for Simultaneous Set of Equations 
i = 1 to 115 

10 


If any |d In ... I >10 calculate -r=rr ; and multiply all D In , 

I |Dln. ..| max 

values by this factor 


19. Test on Condensation Products 
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1 


le 

oe 


Continue otherwise go to step 18 


1 - 


1 - 


1 - 


lU 

ou 

IX 

ox 

lY 

OY 

IZ 

OZ 


DQi 


I < 5 X 10-’ 

< 5 X 10'’ 

< 5 X 10*^ 

< 5 X 10-’ 

< 5 X 10-’ 
5 X 10'^ 


Continue otherwise go to step 18 
Continue otherwise go to step 18 
Continue otherwise go to step 18 
Continue otherwise go to step 18 
Skip to step 19 otherwise go to step 18 


In Pi = InPi + Din Pi 
In A = In A + D In A 
In T = In T + D In T 
Return to step 11 


Determine Pi and test: If Pi = 0 set Pi = 1 x 10 
Determine A 
Determine T 


If 



InPlOl <0 

call for No. Ill and continue 



In Pi < 0 

call for No. 112 and continue 


If 



If 



In P^ < 0 

call for No. 113 and continue 
In Pi < 0 

call for No. 114 and continue 
In Pi < 0 

call for No. 
to step 1 1 


If all values are > 0 go to step 27 


115 and return 



20. Entropy of the Combustion Products [cal/g*K] 


21. Test on Tolerance for Constant Entropy Process at Frozen Equilibrium 
Condition 

22. New Estimates on Temperature for Constant Entropy Process at Frozen 
Equilibrium Condition 


23. Partial Pressures of the Combustion Products for Frozen Equilibrium 
i = 1 to 110 

24. Sum of the Partial Pressures for all Gaseous Combustion Products [atm] 


25. Simultaneous Set of Equations Solving for Derivatives DPP (P = Constant) 

a. Mass Balance Equations 

For definition of DH^ see equation 16a 


122 



no 


R • E 

i = 1 


CNi . In 


CNj • PO 
CNSMCP 


2 

A 


If 



<5x10-’ 


skip to step 23 otherwise continue 


D InT = 


(SO - S) 



In T = In T + D In T 
determine T and return to step 11 


P - CNj-PO 
i “ CNSMCP 


no 

PSUM = Pi 
i= 1 


no 

o = Y. Ci • Pi[ Ci • DPPlOl + Fi • DPP102 + Hi • DPP103 + HEi • DPP104 
i = 1 

+ ANi • DPP105 + 0i • DPP106 + Ui • DPP107 + Xi • DPP108 + Yi • DPP109 

115 

+ Zi- DPPllO + DHi] + Y Ci-Pi-DPPi -,A- IC-DPA 

i= 111 

no 

0 = Y Fi-Pi[Ci- DPPlOl + Fi-DPPlOZ + Hi-DPP103 + HEi-DPP104 
i = 1 

+ ANi- DPP105 + 0i-DPPlO6 + Ui- DPP107 + Xi' DPP108 + Yi*DPP109 

115 

+ Zi . DPPllO + DHj] + Y Pi- Pi- DPPi - A- IF - DPA 

i= 111 
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110 

O = X Hi.Pi[Ci.DPP101 + Fi-DPP102 + Hi-DPP103 + HEi-DPP104 
i= 1 

+ ANi • DPP105 + 0i • DPPI 06 + Ui • DPP107 + Xi • DPP108 + Yi • DPP109 

115 

+ Zi- DPPllO + DHi] + X Hi - Pi DPPi - A • IH • DPA 

i= 111 

110 

O = X HEi • Pi[Ci • DPPlOl + Fi • DPP102 + Hi • DPP103 HEi • DPP104 
i= 1 

+ ANi • DPP105 -1 0i • DPP106 -1- Ui • DPP107 -1 Xi • DPP108 + Yi • DPP109 

115 

Zi • DPPllO + DHi] + X HEi • Pi • DPPi - A • IHE • DPA 

i= 111 

110 

O = ANi -Pi[Ci -DPPlOl + Fi -DPP102 -1 Hi • DPP103 + HEi • DPP104 
i= 1 

+ ANi-DPP105 + 0i-DPPlO6 +Ui-DPP107 + Xi-DPP108 +.Yi-DPP109 

115 

+ Zi • DPPllO + DHi] + X, ANi • Pi • DPPi - A • IN • DPA 

i = 111 

110 

0 = X Oi • Pi[Ci • DPPlOl + Fi • DPP102 + Hi • DPP103 + HEi • DPP104 
i = l 

+ ANi • DPP105 + 0i • DPP106 + Ui • DPP107 + Xi • DPP108 + Yi • DPP109 

115 

+ Zi- DPPllO DHi] + X 0i- Pi -DPPi - A- 10 -DPA 

i = 111 

110 

0=2] Ui -Pi[Ci- DPPlOl + Fi-DPP102 + Hi- DPP103 + HEi- DPP104 
i = 1 

+ ANi • DPP105 + 0i - DPP106 -I- Ui - DPP107 + Xi - DPP108 + Yi • DPP109 

115 

+ Zi- DPPllO + DHi] + X Ui-Pi-DPPi - A-IU-DPA 

i= 111 
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25 a. (Continued) 


Equilibrium Equation for Condensed Products 
i = 111- to 115 


Total Pressure Equation 


126 



110 

= £ Xi-Pi[Ci-DPP101 + Fi-DPPlOZ + Hi-DPPlOS + HEi-DPP104 
i = l 

+ ANi • DPP105 + 0i • DPP106 + Ui • DPP107 + Xi • DPP108 + Yi • DPP109 

115 

+ Zi • DPPllO + DHi] + X Xi • Pi • DPPi - A • IX • DPA 

i= 111 

110 

= Z Yi-Pi[Ci*DPP101 + Fi-DPPlOZ + Hi-DPP103 + HEi-DPP104 
i= 1 

+ ANi • DPP105 + 0i • DPP106 + Ui • DPP107 + Xi • DPP108 + Yi • DPP109 

115 

+ Zi • DPPllO + DHi] + Z Yi • Pi • DPPi - A • lY • DPA 

i = lll 

110 

= y Z. • P.[C. • DPPlOl + F. • DPP102 + H. ' DPP103 + HE. • DPP104 

^ 1 i*- 1 1 1 1 


+ AN- DPP105 + 0J-DPP1O6 + Uj'DPP107 + Xj-DPP108 + Y- DPP109 

115 

+ Zi- DPPllO + DHil + y Zi • Pi DPPi - A • IZ • DPA 

i= 111 

= Ci • DPPlOl + Fi • DPP102 + Hi • DPP103 + HEi • DPP104 + ANi • DPP105 
+ 0i-DPPlO6 -I- Ui- DPP107 -I- Xi • DPP108 -f Yi-DPP109 + Zi-DPPllO-H DH 

no 

= Z PiLCi • DPPlOl -H Fi-DPP102 -f Hi- DPP103 -f HEi- DPP104 
i=l 

+ ANi • DPP105 -^ 0i • DPP106 + Ui • DPP107 -h Xi - DPP108 + • DPP109 

+ Zi - DPPllO + DHi] 


b. Derivatives for the Gaseous Combustion Products 
Equilibrium Equation for Gaseous Compounds 
i = 1 to 100 

For definition of DHj see equation 16a 


26. Simultaneous Set of Equations Solving for Derivatives DTP^ 
(T = Constant) 

a. Mass Balance Equations 
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DPPi = Ci • DPPlOl + Fi • DPP102 + Hi - DPP103 + HEi • DPP104 + ANi • DPP105 
+ 0i ■ DPP106 + Ui • DPP107 + Xj • DPP108 + Yi • DPP109 + Zi- DPPllO 
+ DHi 


115 110 

0 = - ^ Ci-Pi+ X Ci-Pi[Ci-DTP101 + Fi-DTP102 + Hi-DTP103 

i = 1 i = 1 

+ HEi • DTP104 + ANi ' DTP105 + 6i • DTP106 + Ui • DTP107 + Xi • DTP108 

115 

+ Yi • DTP109 + Zi • DTPllO] + Y Ci-Pi-DTPi 

i = 111 

115 110 

0 = - Y Fi-Pi+ Y Fi •Pi[Ci- DTPlOl + Fi -DFP102 + Hi -DTP103 
i = 1 i = 1 

+ HEi • DTP104 + ANi ' DTP105 + Gi • DTP106 + Ui • DTP107 + Xi • DTP108 

115 

+ Yi • DTP109 + Zi • DTPllO] + Y Fi-Pi-DTPi 

i = 111 

115 no 

0 = _ ^ Hi-Pi+ Y Hi •Pi[Ci- DTPlOl + Fi-DTP102 + Hi- DTP103 
i = 1 i = 1 

+ HE. • DTP104 + ANi • DTP105 Gi • DTP106 + Ui • DTP107 + Xi - DTP108 

115 

+ Yi-DTP109 + Zi- DTPllO] + Y Hi-DTPi 

i = 111 

115 no 

0 = - Y HEi- Pi + Y HEi - Pi[Ci- DTPlOl + Fi-DTP102 + Hi- DTP103 
i = 1 i = 1 

+ HE- - DTP104 + ANi - DTP105 + Gi - DTP106 + Ui - DTP107 + Xi • DTP108 

+ Yi - DTP109 + Zi • DTPllO] + Y HEi - DTPi 

i = 111 
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115 no 

0 = - 5] ANi-Pi+ Y ANi-Pi[Ci'DTP101 + Fi-DTP102 + Hi‘DTP103 

i = 1 i = 1 

+ HEi • DTP104 + ANi • DTP105 + 0i • DTP106 + Ui • DTP107 + Xi • DTP108 

115 

+ *Yi • DTP109 + Zi -DTPllO] + Y ANi-DTPi 

i = 111 

115 no 

6 = - Y 6i‘Pi+ Z Oi • Pi[Ci • DTPlOl + Fi « DTP102 + Hi • DTP103 
i = 1 i = 1 

+ HEi • DTP104 + ANi ’ DTP105 + 0i • DTP106 + Ui • DTP107 + Xi • DTP108 

115 

+ Yi-DTP109 + Zi -DTPllO] + Y Oi-DTPi 

i = 111 

115 no 

0 = - 5^ .Ui-Pi+ Y Ui • Pi[Ci -DTPlOl + Fi • DTP102 + Hi • DTP103 

i = 1 i = 1 

+ HEi • DTP104 + ANi - DTP105 + 0i - DTP106 + Ui - DTP107 + Xi - DTP108 

115 

+ Yi - DTP109 + Zi • DTPllO] + Y Ui-DTPi 

i = 111 

115 no 

0 = - Y Xi-Pi + Y Xi • Pi[Ci- DTPlOl + Fi- DTP102 + Hi- DTP103 
i = 1 i = 1 

+ HEi - DTP104 + ANi ‘ DTP105 + 0i • DTP106- + Ui - DTP107 + Xi - DTP108 

115 

+ Yi - DTP109 + Zi - DTPllO] + Y Xi - DTPi 

i = 111 

115 no 

0 = - Y Yi-Pi+ Y Yi -Pi[Ci- DTPlOl + Fi-DTP102 + Hi- DTP103 
i = 1 i = 1 

+ HEi • DTP104 + ANi ' DTP105 + Oi - DTP106 + Ui - DTP107 + Xi - DTP108 

115 

+ Yi- DTP109 + Zi-DTPllO] + Y Yi-DTPi 

i = ,11 1 
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26 a. (Continued) 


Equilibrium Equation for Condensed Products 
i = 111 to 115 


b. Derivatives for the Gaseous Combustion Products 
Equilibrium Equation for Gaseous Compounds 
i = 1 to 100 


27. Separation of Mole Numbers and Partial Pressures Which Will be Used 
for Further Calculations [mole] 
i = 1 to 115 


28. Total Mole Number of the Combustion Product [mole] 


Total Mole Number of the Combustion Product Minus the Condensed 
Phases [mole] 

29. Test on Desired Calculation Procedure During Expansion 


115 no 

’ 0 = - Y, Zi-Pi+ Y Zi • Pi[Ci • DTPlOl + Fi • DTP102 + Hi • DTP103 
i = 1 i = 1 

+ HEi • DTP104 + ANi • DTP105 + 0i • DTP106 + Ui • DTP107 + Xi • DTP108 

115 

+ Yi • DTP109 + Zi • DTPllO] + Y Zi-DTPi 

i = 111 

O = C. • DTPlOl + Fj • DTP102 + • DTP103 + HE- • DTP104 + AN- • DTP105 

+ 0i • DTP106 + Ui • DTP107 + Xi • DTP108 + Yi • DTP109 + Zi • DTPllO 
DTPi = Ci • DTPlOl + Fi • DTP102 + Hi • DTP103 + HEi ' DTP104 + ANi ‘ DTP105 
+ 6i • DTP106 + Ui • DTP107 + Xi • DTP108 + Yi • DTP109 + Zi • DTPllO 


CNi = Pi 


115 

CNS = Y 
i = 1 

no 

CNSMCP = Y CNi 
i = 1 

(a) If shifting equilibrium calculation is desired set SHIFT = 1 and replace 
enthalpy equation in simultaneous set of equations by entropy equation 

(b) If frozen equilibrium calculation is desired set FROZ = 1 

(c) If combination of shifting and frozen equilibrium is desired test 

If T > TL set SHIFT = 1 
If T < TL set FROZ = 1 


If SHIFT si go to step 25 
If FROZ si go to step 31 
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31. For Frozen Equilibrium the Mole Numbers do not Change 


» 


32, Specific Heat at Constant Pressure for Frozen Equilibrium [ cal /mole ®K] 


33. Enthalpy of the Gaseous Combustion Products for Frozen Equilibrium 
[cal/g] 

34. Change of Molecular Weight as a Function of Temperature at Constant 
Pressure for Shifting Equilibrium [ - ] 

35. Change of Molecular Weight as a Fiinction of Pressure at Constant 
Temperature for Shifting Equilibrium [ - ] 

Coefficients of Derivative Compare with Pressure Row in the Matrix 
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Set DPM = 0 
DTM = 0 


CPP = 


115 

Y, CNi-CPi 
i = 1 

no 

I CNi 
i = 1 


r 115 


H = 



CNi* HTi 


A 


DPM = DP A 


and skip to step 37 


dtm 1 

no 

Y, Pi [Ci • DTPlOl + • DTP102 + • DTP103 

i = 1 + HEi • DTP 1 04 + ANi * DTP 1 05 + 0i * DTP 1 06 

+ Ui • DTP107+ Xi • DTP108+ • DTP109 

+ Zi • DTPllO] 
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36. Specific Heat at Constant Pressure for Shifting Equilibrium [cal /mole “K] 
Coefficients compare with enthalpy row in the matrix 


37. Specific Heat at Constant Volume [cal /mole *K] 


38. Specific Heat Ratio at Constant Entropy [ - ] 


39. Mole Ratios of the Combustion Products [ - ] 
i = 1 to 115 


40. Average Molecular Weight Per Mole of Gaseous Products [g/mole] 
The volume of the condensed phases is considered negligible com- 
pared to the gaseous compounds. 
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CPP = 


110 

Z CNi-T 
i = 1 


- R • T 


;r,(s), . 


DPA + R • T 


:? .(f ), - 


+ R • T 

no 


Z (^).- CNi ■ Ci ■ DPPlOl + X ^ (^)i ■ DPP102 

j(^)l CNi ■ HEi ■ DPP104 

{¥^\ CNi • ANi ■ 0PP1O5 + ^X^ (^)i 

(-^). CNi ■ 'Ji • I’PPiO’ * ,‘f ^ {-Ft), ■ 

. 2°^ (i^)l' CNi-yi-DPP109 + ^X^ (^)i' ^ ‘ 


DPP108 


DPPllO 


CVV = CPP - R 


CRS = 


(1 - DPM)^ 
(1 + DTM) 

CPP 

R 


CPP 

R 


[ 1 + DTM] - [ 1 - DPM]' 


AMNRi= 


CNi 


X CNi 

i = 1 


AMW = AMWR 

I CNi 
i = 1 
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41. Convergence 


Delete 
Delete I 1 


(■■tf 


during the first calculation loop 


during the following calculation loops 


42. Specific Gas Constant [cal/g°K] 
R= 1.98726 [cal/mole °K] 


43. Local Speed of Sotind [ft/sec] 
CK3 = 37.42098622 


44, Local Velocity [ft/sec] 

CK4 = 52. 92126623 

1 

— • y HT^ * CNi from simultaneous set for shifting equilibrium 

A . ^ ^ 

1=1 

1 

— • Y, HTi • CNi from step 33 for frozen equilibrium 
^ i = 1 

45. Local Specific Impulse [sec] 

KE = 32. 1740 [ft/sec^l 


46 • Test 

a. For chamber calculation determine the expansion ratio 

b. For throat area calculation set 


c. For all consecutive calculation loops 
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2 


CONVER = \ 1 


no 
Z Pi 

i = 1 
PO 




RS = 


AMW 


VS = CKsn/CRS • RS • T • KE 


V = CK 


:y|io - X . 5] ^ HTi . CNi) 


KE 


AISP = 


V 

KE 


EXPR = 


ACC 

AS 


PS = PO 


TS 


and skip to step 64 
T RSS = RS 


VSS = VS 


CRSS = CRS 


and go to step 47 


go to step 47 
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47. Local Expansion Ratio [ - ] 


48. Weight Flow Rate [lb /sec] 
CKi = 0. 01049466 


49. Characteristic Yelocity [ft/sec] 
CK 2 = 14.696 [Ib/in^atm] 


50. Local Thrust Coefficient [ - ] 


51. Local Thrust [lb] 
CK 2 = 14. 696 


52. Thrust at Altitude [lb] 

CK 2 = 14. 696 

53. Thrust at Vacuum [lb] 

CK 2 = 14. 696 

54. Thrust Coefficient at Altitude [ - ] 
CK 2 = 14. 696 


55. Thrust Coefficient at Vacuum [ - ] 
CK 2 = 14. 696 


56. Real Thrust Coefficient at Altitude [ - ] 


57. Real Thrust Coefficient at Vacuum [ - ] 


58. Real Thrust at Altitude [lb] 
CK, = 14.696 


59 . Real Thrust at Vacuum [ lb] 
CK 2 = 14.696 
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_ RS • T • VSS • PS 
" V • PO • RSS • TS 


w = 


_ PO • V • AS • EXPR • CKi 
RS • T 


PC-AS-KE-CK, 

c-= w ' 


CF = 


AISP • KE 
C- 


FT = AS • PC • CF - CK2 


FTA = FT - EXPR • AS • (PAM - PO) • CK2 


FTV = FT EXPR • AS • PO • CK2 


CFA = 


FTA 

AS • PC • CK2 


CFV 


FTV 

AS • PC • CK2 


CFRA = ETACFV • CFV 


EXPR- 


PAM 

PC 


CFRV = ETACFV- CFV 
FTRA = CFRA • AS - PC • CK2 


FTRV = CFRV ■ AS - PC • CK2 



60. Real Weight Flow Rate [lb /sec] 

61. Real Characteristic Velocity [ft/sec] 

62. Real Specific Impulse at Altitude [sec] 

63. Real Specific Impulse at Vacuum [sec] 

64. a. Calculate X value for plotting purpose 

b. Store: 

Options 

65. Test: If lA = 0 conditions in the throat area are calculated 

If LA 1 conditions except the throat area are calculated 

66. Test: If IB = 0 print the results of the previous calculation and 

determine the approximate throat parameters 


If IB = 1 
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WR = ETAW • W 


C^<R = ETAC>!< • C* 


AISPRA = 


FTRA 

WR 


AISPRV = 


FTRV 

WR 


If V < VS calculate X = EXPR - 1 
If V > VS calculate X = -(EXPR - 1) 

During the first calculation loop set: TC = T 

SO = S 

During every calculation loop store the results as indicated on page 107 
for printing purpose 


If LA = 0 skip to step 66 
If LA = 1 go to step 69 

If IB = 0 PRINT, then calculate 

/ CRS 

r 2 • * 

[ crIttJ 

^ = "'^[crItt] 

Set IB = 1 and return to step 11 
If IB = 1 go to step 67 
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67. Test on Local Velocity and Speed of Sound 


68. Test: IC is an Input Parameter Indicating 

a. IC = 0 only chamber, throat and exit conditions are considered 

b. IC = 1 intermediate locations between chamber, throat and exit 

area are considered 

69. Test: Determination of whether the point of calculation is upstream 

or downstream of the throat area 
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If 


1 - 


vs 


PO = PO 


> 0. 2 X 10"* calculate 


R- T 


2 • AMW 
and return to step 11 
V 


(CRS + 1) 


If 


1 - 


VS 

Set: IB = 0 


< 0. 2 X 10-'* 


lA = 1 
PRINT 

and go to step 68 


If IC = 0 calculate 

oo _ „-[ (0. 867- CRS + 0. 272)lnEXPRL -1- (0.699- CRS + 0.442)] 

Ir vJ — * hj 

and return to step 11 

If IC = 1 calculate 

PO = PC - 10 

and return to step 11 

If IC = 0 test 

If VS > V skip to equation PO in step 70 
If VS < V go to step 70 
If IC = 1 skip to step 71 
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70. Test on Required Expansion Ratio 


71. Test on Expansion Ratio 

The test on VS and V shall avoid an iteration on an expansion ratio in 
the subsonic range 


72. 


Test to Define Assumed Static Pressure Differences for Various 
Ranges during the Expansion Process 


73. Static Pressure in the Nozzle Between Chamber Pressure and 1 atm 
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If 


1 - 


EXPR 

EXPRL 


> 1 X 10“^ calculate 


PO = PO 


HO 


H + 


PS^ • CRSS • RS^ • 


1 + 


RS 


•T^l - 


RSS • TS • 2 • PQ2 • EXPRL2, 
CRSS-PS^ - RS -T 


) 


CRS • PO^ • RSS • TS • EXPRL^, 




and return to step 1 1 


If 1 - 


EXPR 


EXPRL I ^ 
set lA = 0 


< 1 X 10 


-3 


PRINT 

and go to new input reading 


If EXPR > EXPRL 
Test: If vs’ > V 
PRIHT 

and go to step 72 
If VS < V continue 

PS^ • CRSS • RS^ ' T^ 

RSS • TS • 2 • EXPRL2 • PO^ 
CRSS • PS^ • RS • T 
RS • PQ2 . RSS • TS • EXPRL^j 

set IC = 0 

and return to step 1 1 
If EXPR < EXPRL 
PRINT 

and go to step 72. a. 

(a) If ID = 1 skip to step 79 otherwise continue 

(b) If IE = 1 skip to step 77 otherwise continue 

(c) If IF = 1 skip to step 75 otherwise go to step 73 



PO = PO 


1 + 


■(. 


HO - I H + 


RS 


•T^l - 


PO = PO - lO 
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74. Test 


75, Static Pressure in the Nozzle Between 1 atm and 0. 2 atm 

76. Test 


77, Static Pressure in the Nozzle Between 0, 2 atm and 0. 1 atm 

78. Test 


79, Static Pressure Between 0. 1 atm and 0,02 atm 

80, Test 


81. Test: If IG 
If IG 


0 maximum specific impulse can be calculated 

1 desired thrust level for a fixed mixture ratio is 
calculated 
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If PO > 1 
return to step 1 1 
If PO < 1 
set PO = 1 
IF = 1 

return to step 11 
PO = PO - 0. 2 


If PO = 0. 2 

set IE = 1 PO = 0. 3 
IF s 0 

return to step 11 

PO = PO - 0. 1 

If PO =: 0. 1 
set ID = 1 
IE = 0 

return to step 11 

PO = PO - 0. 02 

If PO < 0. 02 
set-IA = 0 
ID = 0 

go to new input reading 

If IG = 0 go to step 82 
If IG = 1 skip to step 103 



82. Test: If IH = 0 no maximum specific impulse is considered 

If IH = 1 maximum specific impulse is calculated for a constant 
chamber pressure and varying mixture ratio (exit 
pressure equals ambient pressure) 

83. Test: During two calculation loops II = 0, afterwards II = 1 

84. New Mixture Ratio for the Approach of Maximum Specific Impulse 
AMRWD is Loaded, but an Assumed Value 

85. 

86. Determination of the Slope for a Tangent through Two Consecutive 
Points on the Specific Impulse Versus Mixture Ratio Curve 

87. Test: Calculation for maximum specific impulse 

Check whether maximum has been obtained 

88. Test: Slope of tangent is not small enough K = 0. 001 

Slope of tangent is within required tolerance K = 0. 001 

89. New Mixture Ratio for the Approach Towards Maximum Specific 
Impulse 
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If IH = 0 go to new input reading 
U IH = 1 set; AISPZ = AISPY 
AISPY = AISPX 
AISPX a AISP 


AMRWZ a AMRWY 
AMRWY a AMRWX 
AMRWX a AMRW 


and go to step 83 


If II = 0 go to step 84 
If II = 1 skip to step 86 

AMRW = AMRW + AMRWD 


Set PO = PC 

II a 1 after second calcvilation loop and return to step 3 


tan a a 


AISPX - AISPY 
AMRWX - AMRWY 


If IJ a 0 go to step 88 
If U a 1 skip to step 93 


If tan a < |k| skip to step 90 
If tan a > |K| go to step 89 


AMRW a 


AMRWY + AMRWX 
2 


AISPY - AISPX 
AMRWY - AMRWX 


AISPZ - AISPY 


AISPY - AISPX 


AMRWZ - AMRWY ~ AMRWY - AMRWX 
AMRWZ - AMRWX 


set PO a PC 

and return to step 3 
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90. Store Maximum Values for Constant Chamber Pressure 


91 . New Mixture Ratio to Check on the Obtained Maximum Specific Impulse 
Value K = 0.001 


92. 

93 . Test: Maximum Specific Impulse is Obtained 

Maximum Specific Impulse is not Within the Required Tolerance 


94 . Repeat Iteration Using the Mixture Ratio Calculated in Step 91 

95 . Print the Maximum Values Stored in Step 90 

96 . Initialize 

97 . Test: Determine the derivation of the loaded mixture ratio from the 

one for maximum specific impulse during the first calculation 
loop 

Consecutive Passes 

98. Calculate the Deviation of the Original Loaded Mixture Ratio from 
the One for Maximum Specific Impulse at Loaded Chamber Pressure. 
Sign is Important for Step 101 
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AISPM = AISP 
AMRWM = AMRW 
PCM = PC 

AMRW = AMRWM + • K 

[tan or| 

Set IJ = 1 PO = PC and return to step 3 

If tan a changes sign and 
If AISPM > AISP go to step 95 
If tan a does not change signs or 
U AISPM < AISP go to step 94 

Set IJ = 0 PO = PC and return to step 89 

PRINT [maximum values for constant chamber pressure] 

Set II = 0 
IJ = 0 

If IK = 0 go to step 98 


If IK = 1 skip to step 100 


PERCT = 


AMRWM - AMRWL 
AMRWM 
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99. Set 


100. Test: If IL = 0 no further calculation for desired thrust level is > 

considered 

If IL s 1 calculation for a desired thrust level deviating in 

mixture ratio by PERCT from the one for maximum 
specific impulse 

101. Calculate Mixture Ratio Which Deviates by PERCT from the One 
for Maximum Specific Impulse at Presently Considered Chamber 
Pressure 

102. Set Switch to Avoid Maximum Specific Impulse Calculation for the 
Following Calculation Pass 


/ 

103. Test: Desired thrust level has been obtained 

Ks 0. 001 


Desired thrust level has not been obtained 


104. Test: First calctilation loop 

Consecutive passes 

105. Set 
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If IL = 0 go to new input reading 
If IL = 1 go to step 101 


AMRW = AMRWM - AMRWM- PERCT 


Set IG = 1 
PO = PC 

and return to step 3 


If 


FTRA 

FL 


< K 


set IG = 0 


IK = 0 


IM = 0 
PRINT 


and go to new input reading 


If 


FTRA 

FL 


> K 


go to step 104 


If IM = 0 go to step 105 
If IM = 1 skip to step 108 


Set: PX = PC 

FX = FTRA 



106. Calculate New Combustion Chamber End Stagnation Pressure 
(First Calculation Loop) 

107. Set 


108. Set 

109. Calculate New Combustion Chamber End Stagnation Pressure 
(Consecutive Passes) 

110. Set 


111, Test: If maximum specific impulse values are considered set IG 
If constant mixture ratio is considered set IG = 1 


PC = PC 


PO = PC 
IM = 1 
and skip 

PY =PX 
PX = PC 

PC = PX 


PO = PC 

If IH = 1 
If IH = 0 


FL 

FTRA 


to step 111 


FY = FX 
FX = FTRA 


PX - PY 
FX - FY 


(FL - FX) 


set IG = 0 
set IG = 1 


and return to step 3 
and return to step 3 





APPENDIX H 


ENGINEERING FLOW CHART 
FOR THE 

CALCULATION PROGRAM 


The symbol in each box refers to the 
proper equation or test in Appendix G 


























I 

I 


F GENERAL COMBUSTION PROGRAM 
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